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ABSTRACT 

We present the photometric catalogs for the star-forming cluster NGC 602 
in the wing of the Small Magellanic Cloud covering a range of wavelengths from 
optical {HST/ACS F814W and SMARTS/ANDICAM F, I) to infrared 

{Spitzer /IRAC 3.6, 4.5, 5.8, and 8 /xm and MIPS 24 /xm). Combining this 
with IRSF (InfraRed Survey Facility) near-infrared photometry (J, H, Kg), we 
compare the young main sequence (MS) and pre-main sequence (PMS) popu- 
lations prominent in the optical with the current young stellar object (YSO) 
populations revealed by the infrared (IR). We analyze the MS and PMS popu- 
lation with isochrones in color-magnitude diagrams to derive ages and masses. 
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The optical data reveal ~ 565 PMS candidates, low mass Stage ill YSOs. We 
characterize ~ 40 Y SOs by fitting their sp ectral energy distributions (SEDs) 



to a grid of models ( iRobitaille et al. 1120071 ) to derive luminosities, masses and 



evolutionary phase (Stage l-lll). The higher resolution HST images reveal that 
~ 70% of the YSO candidates are either multiples or protoclusters. For YSOs 
and PMS sources found in common, we find a consistency in the masses derived. 
We use the YSO mass function to derive a present-day star-formation rate of 
~ 0.2 — 1.0 Mq yr^^ kpc~^, similar to the rate derived from the optical star 
formation history suggesting a constant star formation rate for this region. We 
demonstrate a progression of star formation from the optical star cluster center 
to the edge of the star forming dust cloud. We derive lifetimes of a few 10^ years 
for the YSO Stages I and ll. 

Please note that color images are compressed for space. Please contact the 
first author for superior, high resolution versions. 

Subject headings: galaxies: star cluster — Magellanic Clouds — open cluster and 
association: NGC 602 — stars: formation, pre-main sequence — ISM: individual 
(N90) 



1. INTRODUCTION 



The Small Magellanic Cloud (SMC) is a valuable astrophysical laboratory for under 
standing the processes of star formation in a galaxy that is extremely dif ferent from the 
Milky Way. In particu lar, it has a subsolar chemical abundance of Z~ 0.004 (IRoUeston et al. 



1999 



Lee et all 20051) and a low dust-to-ga s ratio of ~ 1/30 Milky Way in the diffuse in 



terstellar medium ( IStanimirovic et al.ll2000l ) and ~ 1/6 in star-forming regions (iBot et al. 



20071 ). The SMC's lack of organized rotation means that star format ion is predominantly 
driven by a combinatio n of tidally-induced cl oud-cloud interactioiis (|Zaritskv et al.l l2000l : 
Zaritsky fc Harris! |2004[) and shell forma tion ( iHatzidimitriou et al.l |2005[ ) . With its close 
proximity (60.6 kpc; iHilditch et al.ll2005l ). the SMC is uniquely suited to detailed investi- 
gation of the stellar content, down to the sub-solar mass regime, in the youngest and most 
compact star clusters. 

The small, yo ung star cluster NGC 602, associated with the highly structured H ll region 
N90 (jHenizd 119561 ). is one of the most interesting star forming regions in the SMC. Located 
at the boundary between the SMC wing and the Magellanic Bridge, at the intersection of 
three Hi shells, it is probably the result of the interaction o f two expariding s hells of Hi 
that occurred approximately ~ 7Myr ago ( iNigra et al.ll2008l ). iNigra et al.l (120081 ) also show 
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that the N90 region is quiescent, with neghgible Ha shell expansion velocities. Propagation 
of star formation is most likely driven by radiation with stars forming along the edges of 
the photodissociation region. Star formation started approximately 4 Myr ago, with the 
formation of the central cluster, and gradua ll y propagated towards t he outskirts, where 
star formation continues (jCarlson et al.l 120071 ). iGouliermis et al.l ( 120071 ) provide a list of 22 
candidate Young Stellar Objects (YSOs) in the outskir ts of N90. The optical cluster's Mass 
Function (MF) is con si stent with a standard Salpeter ( 1Salpeteiill955l ) Initial Mass Function 
(jSchmalzl et al.ll2008l ). ICignoni et al.l ( l2009l ) reconstruct a complete star formation history 
for the optical population and find that the pre-main sequence (PMS) is not more than 
~5 Myr, and the star formation rate (SFR) has reached a maximum in the last 2.5 Myr. 

We employ a new panchromatic approach to extend the star formation history analysis 
to the present day. Population identification via Color-Magnitude Diagram (CMD) analysis 
informs us of cluster-scale star formation, with optical revealing a bright main sequence (MS) 
and a young PMS and infrared (IR) CMD highlighting the youngest embedded sources. Using 
spectral energy distribution (SED) analysis, for which we combine Hubble Space Telescope 
(HST) optical to probe central stellar sources and Spitzer Space Telescope (Spitzer) IR to 
probe circumstellar disks and envelopes, we characterize star formation in the NGC 602 
region (Figure [1]) on the scale of single sources and protoclusters. 

The structure of this paper is as follows. We describe observations and data reduction 
in Section [21 In Section [31 we discuss what can be gleaned from the IR and optical data 
separately via CMD analysis. We describe the application of the SED fitter in Section [H 
and present the results of our SED analysis of 77 sources, combining optical to mid-IR data, 
for YSOs (|Ml KM KB and SZj) and other types of IR sources flilSj) . We address the YSO 
Mass Function and Star Formation Rate (SFR) in Section [51 and the spatial and temporal 
distribution of young sources in Section [HI 



2. OBSERVATIONS AND DATA REDUCTION 

Observations of NGC 602 (R.A. = 01^29'^3P, Dec. = -73°33'15", J2000) were taken 
from the optical to the mid- infrared (Figure [JJ. Optical photometry (Section 12. ip reveals 
the naked MS and PMS populations, and we determine their ages through CMD analysis. 
Accompanying HST Advanced Camera for Surveys (ACS) imaging offers a high- resolution 
view of both the nebular structure and the physical distribution of the stars. Spitzer Infrared 
Array Camera (IRAC) data (Section 12. 2p reveal the YSOs. Relying completely on Spitzer 
photometry and imaging, however, would introduce significant contamination from back- 
ground galaxies, contamination which can be avoided through comparison to ACS images 
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(See Section 14.8.21 and iNigra et al.l l2008l ). Likewise, optical data alone miss the youngest 
sources. We characterize YSOs and other infrared sources by combining optical, near in- 
frared (NIR), and IR photometry (Section 12.31) and fitting the resultant SEDs to those of 
numerically modeled YSOs, catalog naked stars, and galaxies. These sources are referenced 
below with prefixes Y (YSO candidate), G (galaxy), K (naked star), S (stellar in images but 
not fit), and U (unclassified). 



2.1. Optical 

We use optical data primarily from the HST/ACS Wide Field Channel {WFC) in 
F555W (V), F658N (Ha), and F814W (i). Some bright Red Giant and upper-MS stars 
are saturated in our ^CS" images. We correct for this by adding ground-based observations 
from the SMARTS (Small and Moderate Aperture Research Telescope System) ANDICAM 
instrument on a 1.3 m telescope in Johnson- Cousins Fand / bands. The log-book of optical 
observations is pres ented online in Table [U The final optical catalog (Table HI also used in 



Cignoni et al.ll2009l ) includes 4535 stars, from a combination of these HST and ground-based 



data. 



2.1.1. HST Data 

Five long exposures were taken through filters F555W (~ V) and F814W (~/) using 
HST/ACS/WFC (JaUe^ i/^T Proposal 10248, RI. Nota). The dither pattern allowed hot 
pixel removal and filled the gap between the two 2048 x 4096 pixel detectors. This dithering 
technique also improved both Point Spread Function (PSF) sampling and photometric accu- 
racy by averaging flat-field errors, and by smoothing over the spatial variations in detector 
response. The entire data set was processed with the CALACS calibration pipeline, and the 
long and short exposures for each filter were separately co-added with the MULTIDRIZZLE 
package. The final exposure times of the deep combined F555W and F814W are 2150 s and 
2265 s respectively. The images cover a region of 200" x 200" corresponding to linear size of 
58 pc X 58 pc at the SMC distance. 

The photometric reduction of all optical images ha s been perforraed wi thin the IRAF0 



environment via the same method used for NGC 330 (jSirianni et al.l |2002| ) and NGC 346 



^IRAF is distributed by the National Optical Observatories, which are operated by AURA Inc., under 
cooperative agreement with the National Science Foundation. 
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( ISabbi et al.l 120071 ). We use P SF fitting and aperture pliotometry routines provided within 
IRAF's DAOPHOT package (IStetsonI 119871 ) to derive accurate photometry of all stars in 
the field. The photometr y is calibrated into the HST Vegamag photometric system, fol- 
lowing the prescription of ISirianni et al.l ( 120051 ). Charg e Transfer Ef f icienc y (CTE) correc- 
tions are included, following the procedure described in ISabbi et al.l ( 120071 ). We consider a 
source photometric in F555W (in F814W) only if its detection is at least Aa (3. Sex) above 
the background level, and its DAOPHOT sharpness parameter is —0.5 < sharpness < 0.5 
(—0.3 < sharpness < 0.3). This allows us to maximize inclusion of real point sources while 
avoiding inclusion of "fuzzy" background galaxies. Our final catalog from WFC includes 
4496 sources. For further sp ecifics of our photora etry, including completeness testing and 
detailed selection criteria, see lCignoni et al.l ( 120091 ). 



2.1.2. SMARTS Ground-Based Data 



A few of bright stellar sources in NGC 602 are saturated in the ^4 CS" images. Photometric 
observations were taken the night of 14 October 2006 (JD: 2454023.67) and include three 
45 second F-band exposures and three 30 second /-band exposure^. Some dithering was 
applied to cover bad pixels. 

We employ IRAF packages XREGISTER, IMSHIFT, and IMCOMBINE to properly 
align and combine the three images in each band. We perform aperture photometry on 
a subset of stars across the FOV, finding the magnitude of each star at aperture radii 
3, 4, ... , 19, and 20 pixels and select an aperture size of 13 pixels (4'.'8) for use in all fur- 
ther photometry. This aperture captures the most stellar flux with the least background 
contamination. We use DAOFIND to construct a photometry list with DAOPHOT pa- 
rameters 0.1 < sharpness < 0.9 and —0.8 < roundness < 0.8, thus eliminating extended 
sources, bad pixels, and cosmic rays from our list. We complete PSF photometry using 
packages PHOT, PSF, and ALLSTAR, requiring a values of 2.3 (in V) and 1.6 (in i). 
We use SMARTS standard-star observations, air mass, and air mass coefficients from the 
same night to calibrate our cluster photometry into the Landolt Vegamag system. See the 
SMARTS/ANDICAM website! for more information. 

Measured magnitud es are then converted from the Landolt to the HST Vegamag system 
using the equations from lSirianni et al.l (120051 ) and used to fill in for sources saturated in our 



imaging and photometry were also performed but are not presented here. Please contact authors 
for further information. 

http://www.astronomy.ohio-state.edu/ANDIC AM 
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HST observations. These are optically bright sources; none are YSO candidates. We replace 
HST data photometry with SMARTS for mp^^^w < 18 and for mpsi4w < 17.5. In Table |2l 
15 sources have SMARTS photometry only and twenty-six have a combination of SMARTS 
and HST photometry. 



2.2. IRAC Mid-Infrared Data 



Five mid-IR images we re taken through e ach of the four Spitzer/IKAC bands (3.6, 4.5, 
5.8, and 8 fim) described in iFazio et al.l ( 120041 ). The data were taken in 30 s High Dynamic 
range mode, using five dithers in the medium cycling pattern (Program ID 125, P.I. Fazio, 
AOR 12485120). For each band the total integration time is 134 s. The Basic Calibrated 
Data (BCD; version Sll.4.0) were used. 

For each channel the IRAC mosaics have been constructed us ing the Spitzer Mosaicker 
software, along with the IRACproc package ( ISchuster et al.ll2006l ). The overlap correction 
module was used to minimize the instrumental offsets between the frames. The total field 
of view (FOV) is approximately 5' x 5', corresponding to a linear size of 88 pcx88pc at the 
distance of the SMC. The native IRAC pixel scale is 1'.'2; through drizzling dithered images, 
we reach a pixel scale of O'.'4/pixel. 

We use IRAF tasks DAOFIND and PHOT to locate and perform photometry on sources 
in each of the mosaics. A Aa cutoff was used in source identification in DAOFIND. Photom- 
etry was performed with PHOT, using a I'.'G (4 drizzled pixels) aperture radius. This is close 
to the IRAC PSF at [3.6] and encircles most of the flux without including nearby sources in 
crowded fields. Background levels are determined and subtraction performed based on the 
background levels in the annuli with inner radius 14.1 pixels and o uter radius 28.2 pixe ls (or 
5'.'64 and 11'.'28), chosen through image examination. Zero-points ( iReach et al. Il2005[ ) were 
applied to normalize photometry to the nominal 12" (30 drizzled pixels, 10 IRAC native pix- 
els) radius and to calibrate to the Vega magnitu d e syst em. For SED fitting (See Section! 



fluxes are calculated according to iReach et al. I ( 120051 ) 



Our deep IRAC photometric catalog (Table |3]) includes 497 sources, each detected in 
at least one IRAC band. Of these forty-eight are photometric in all four IRAC bands and 
forty-one in three bands, using aperture photometry. The other 408 sou r ces ar e photometric 
in only one (164 sources) or two (244 sources) bands. iGouliermis et al.l ( 120071 ) use the same 
IRAC data set and apply PSF photometry to produce a list of 103 sources detected in at 
least three IRAC bands. In our online p hotometric catalog, we include cross identification 
with the 22 YSO candidates reported by IGouliermis et al.l (120071 ). 
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For the purposes of image examination, we supplement our deeper IRAC data with 
photometry from 12 s exposures from the Spitzer Legacy Progr am SAGE-SMC (Su rveying 
the Agents of Galaxy Evolution in the Small Magellanic Cloud; iGordon et al.ll2010l ). While 
this does not significantly increase the depth, the mosaicing of images taken with differing 
camera orientations allows us to create cleaner images with rounder point spread functions. 
We also use SAGE-SMC IRAC photometry for 3 sources within our 3' radius but outside of 
our deeper IRAC FOV. 



2.3. Supplemental Photometry for SED Construction 

We characterize individual sources by constructing SEDs and comparing them to the- 
oretical YSO model SEDs a nd SEDs of catalog naked stars and galaxies (See Section I4.1j 



Robitaille et al. Il2006l 120071 ). We need to include as many photometric data points as pos- 
sible to best constrain fits. For this analysis, therefore, we require sources to have at least 
4 bands of photometry, including [8.0] and at least 2 other IRAC bands. We examine the 
77 IRAC "fitter sources" meeting these criteria. We add optical and NIR data where avail- 
able, but differences between instrument resolution, particularly at extragalactic distances, 
complicate this. Of our 77 fitter sources, 61 lie within the ACS FOV, 6 with no apparent 
optical (point-source) counterpart, 29 corresponding to a single optical source, and 26 with 
multiple optical counterparts. We must handle our optical fluxes carefully. 

Further complication results from the fact that the SED fitter allows input only in spe- 
cific photometric ban ds, including Landolt Fand / and Two-Micron All Sky Survey (2MASS; 



Skrutskie et al.l 120061 ) J, H, and K^. We therefore convert magnitudes from similar native 
bands, introducing an additional element of uncertainty, as conversions are generally cali- 
brated for MS stars rather than PMS or YSO sources. Our treatment of these uncertainties 
is discussed below. Our applied error estimates are appropriate with fluxes tracing contin- 
uous curves and no discontinuities between fluxes from different instruments except where 
fluxes are used as upper limits for the physical reasons explained below. 



2.3.1. Optical 

Comparing deep, high- resolution ACS optical and IRAC images and photometry re- 
quires special care. Several physical sources at the distance of the SMC may fall within the 
1'.'6 IRAC [3.. 6] aperture. Something that appears to be a single source in the IR images, 
may be resolved into ten or more sources in ACS images, with no obvious way to tell which 



- 10 - 



one is or which ones are the primary IR contributer(s). Likewise, faint or embedded sources 
may have no resolvable optical point source within the IR aperture. 

We consider four possible types of IR-optical matches: 

1. Single Optical Point Sources: These are unsaturated in the {HST) optical and have 
only one obvious optical point source. 

2. No Optical: There is no apparent optical counterpart to the IR source. These often 
correspond to dusty structures remmeniscent of the famous pillars and mountains of 



creation (e.g., Thompson. Smith, fc Hesterl 120021 ). 



3. Multiple Optical Sources: Some of these are probable protoclusters, while others are 
superpositions. These may be faint or bright. 

4. Saturated Sources: Optically saturated sources mostly seem to be bright Main Sequence 
stars. In some cases, faint neighbors with poor photometry (often because of diffraction 
spikes) are distinguishable in optical images. 

We want to include light from the same sources in each of the photometric bands to 
the greatest extent possible to produce the most physically meaningful SEDs (see Section 
HJ. We perform aperture photometry with DAOPHOT's PHOT task at the exact IRAC 
coordinates of the 43 fitter sources which lie within our ACS FOV and are not saturated in 
the optical and use a 32 ACS pixel aperture radius to measure the integrated optical flux 
over the 1'.'6 IRAC aperture. For 13 saturated sources, we use the SMARTS ground-based 
photometry (aperture radius 4'.'8) and calculate the flux as with other sources and apply a 
40% uncertainty. One optically saturated source (S456) has no good optical photometry. 

The SED fitter requires the input of fluxes in specific photometric bands, including 
Landolt V and / but not ACS F555W or F814W. We must convert from the ACS to the 
Landolt system and calculate fluxes from magnitudes. To this end, photometry is performed 
in the ACS OBMAG system (PHOT's ABMAG with no zero poi nt). We then c onver t 



to the Landolt photometric system according to the description of ISirianni et al.l ( l2005l ) . 



We calculate the fluxes in mJy, us ing the zero po int fluxes /oy = 3.789 x 10 mJy and 



foi = 2.635 X 10^ mJy adapted from I Crawford! ( 1l975l ). The conversions from ACS to Landolt 



bands are not designed for dusty, young, or low-metalicity sources. We therefore calculate the 
errors in flux by squaring the errors described in the conversion. Optical fluxes are treated as 
upper limits for several YSO candidates, if their high resolution optical morphology indicates 
the presence of multiple optical or zero optical counterparts. Source Y217 has indefinite 
OBMAG in F555W but a measurable magnitude in F814W. We calculate the flux in F814W 
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using the equation specified in ISirianni et al.l ( 120051 ) and assume that the real / flux is not 



more than a factor of two greater. We thus double calculated flux as an upper limit in the 
SED fitter input (Sections 14. II and I 



In multiple optical sources , it is unlikely that all optical sources are significant IR emit- 



ters (i.e. 



Bernasconi fc Maedeiill996l . Lyso oc Myso)- is instead probable that one, two, 
or even zero of the optical sources is "the" infrared emitter. These sources should be treated 
with extra care. The approximate number of optical sources corresponding to a single IR 
source is determined visually, as some optical sources are non-photometric (See Section I^TT]) . 
Optical fiuxes are treated as upper limits for several YSO candidates, based on their high 
resolution optical morphology. For sources that have no optical or are saturated in the op- 
tical but have 4 IR bands to constrain the SED fit, we set the calculated Fand / fiuxes as 
upper limits rather than assigning errors bars. They are noted in Table which provides 
details of the fitter results as well as the number of optical sources. Optical matches to IRAC 
sources are also noted in online photometry Tables |2] and [3l 



2. 3. 2. IRSF Near Infrared Data 



W e use Near-IR photometry from the Magellanic Clouds Point Source Catalog of lKato et al 
( 120071 ). This imaging survey was performed with the Simultanious three-color InfraRed Im- 



ager for Unbiased Surveys (SIRIUS) aboard the InfraRed Survey Facility (IRSF) 1.4 m 
telescope at the South African Astronomical Observatory's Southland. Photometry cov- 
ers the J (Ac = 1.25 /im), H (1.63 /xm), and Ks (2.4 /im) bands. We choose to use this 
IRSF/SIRIUS data for its depth (< 16 mag in Kg compared to < 14 mag in 2 MASS) and 



high s patial resolution {FWHM ~ I'.'l - 1'.'3, compared to 2MASS's 2" - 3"). iKato et al. 



( 120071 ) performed PSF photometry for aperture radius 1".35 to construct their catalog. 



Sixty-four fitter sources have IRSF matches. For use in the SED fitt er, we convert IRSF 
magnitudes to 2MASS fiuxes. We apply magnitude conversions given in iKato et al.l ( 120071 ). 
which require at least 2 bands of IRSF data, then apply 2MASS zero -points of 1.594, 1.024, 
and 0.6668 mJy in J, H, and Ks ( ICohen. Wheaton. fc MegeathI 120031 ) respectively to obtain 
2MASS fiuxes. Unfortunately, uncertainty is always introduced when converting between 
photometric bands, particularly for non-photospheric objects. We estimate a nominal 20% 
uncertainty in fiux. Eight of the 64 matched sources have only one band of IRSF photometry, 
making the usual conversion to 2MASS magnitudes impossible. We use the IRSF magnitude 
as though it were a 2MASS magnitude and account for extra uncertainty introduced by lack 
of color information by increasing our applied error bars from 20% to 40% in fiux. Another 
13 of the 64 matched fitter sources have multiple matches within the IRSF catalog. For 



these, we take the brightest IRSF measurement in each band and treat the flux as an upper 
limit in our SED fitting. 



2.3.3. MIPS Data 



Multiwavelength Imaging Photometer for Spitzer (MIPS; 



Rieke at al. 



2004 



24 fim 



data are taken from the first epoch observations of SAGE-SMC (iGordon et al.ll2010l ). These 
observations were performed in October 2 007. Post-processing is performed as it was for 



the SAGE-LMC MIPS 24 fim data (as in iMeixner et al.ll2006[). Ph otometrv is performed 



using the shape-fitting algorithm STARFINDER (jPiolaiti et al.ll2000[ l. Sources for which we 
include 24 /xm fiuxes have been selected through automated matching to our known IRAC 
sources and verified through careful visual inspection of the 24 /im image. We match eleven 
24 /im sources to sources fit as YSO candidates and seven more to other types of sources 
(See Section H]). 



3. COLOR MAGNITUDE ANALYSIS - STELLAR AND PROTO-STELLAR 

POPULATIONS 

We examine optical and IR data independently before combining them. This strategy 
results in a complete picture of stellar and proto-stellar populations in the region with the 
more evolved populations revealed by optical alone and the least evolved sources appearing 
only in the IR. 



3.1. The Optical CMD 



Examination of the optical CMD (Figure |3]) reveals a quantitative picture of the optical 
populations. In NGC 602, optical point sources comprise MS, Red Giant Branch (RGB) 
evolved, and PMS stars. Old (lower) MS and th eir RGB counterparts are found to be 
primarily background sources (jCignoni et al.l l2009l ). More important for this study, the 
optical CMD reveals young (< 5 Myr) MS and PMS stars. 

Our optical catalog of over 4500 sources is a combination of HST/ACS and SMARTS 



data (Table The applied 



for NGC 602 by Cignoni et al. 



MS ev olutionary tracks (and isochrones) were calculated 



(cf. IChieffil Il989l : iDegFInnocenti et al. 
M)o = 18.9 and reddening E{B 



(m 



(12009 ) work ing from FRANEC stellar evolutionary code 
20081 ) and using Z = 0.004 and distance modulus 
V) ~ 0.08. Padua stellar evolutionary tracks 
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( iFagotto et al.l Il994[ ) are used for MS and RGB populations. The RGB and lower-MS 
populations are best fit with Z = 0.001, consistent with the Magellanic Bridge/Tail field 
population ( IZaritsky &: Harris! 12004 ). supporting the supposition that these are background 
sources. 



We identify PMS sta rs via the color- magnitude selection (yellow line in Figure 3(a)) 
from ICarlson et al.l ( 120071 ) : 



(28 - 4.3 X [mF555W - "^FSUw]) < "^F555W < (5.5 X [mF555W - "^F814W] + 17) (1) 

In Figure [31 r nost of the PMS po pulation lies between the example 1 Myr and 10 Myr 



isochrones. In ICignoni et al.l ( 120091 ) . it is shown that most PMS stars in NGC 602 are less 
than 5 Myr in age, and the star formation rate has been increasing for the past few Myr. 
The star formation rate was ~ 1.5 x 10~^ Mq yr^-"^ from 5 to 2.5 Myr ago and has increased 
to ~ 3 X 10~^ Mq yr~^ in the last 2.5 Myr. Further detailed analys i s of th e star formation 
history of the optical population can also be found in ICignoni et al.l (120091 ). 



3.2. The Infrared Population 

Figure m shows both a CMD and a Color-Color Diagram (CCD) in IRAC featuring the 
77 sources within 3' of cluster center for which we perform detailed SED analysis below. 
SAGE-SMC IRAC sources for the larger 18' x 18' FOV (317 pc x 317 pc; Figure ^ are 
plotted in grey. This background population is quite sparse and is concentrated near the 
zero-color MS. 

Infrared sources in the NGC 602 region may be classified into 3 groups: YSOs, other 
stars, and background galaxies. Significant features are the colorless MS and the primary 



YSO color-magnitude space, centered around [3.6] — [8.0] = 3 (Figure 4(a)). In the color- 
range 0.3 < [3.6] — [8.0] < 3, only 6 sources are brighter than about [3.6] = 13.5, indicating 
that there are virtually no IR-bright evolved sources in the entire region. Six fit stellar 
sources (K181, S213, K225, K232, K441, and K444) have colors 0.6 < [3.6] - [8.0], redward 
of the MS, and are discussed further in Section [4.8.11 Five of the 7 identified galaxies on the 
CMD lie in the color range 2.5 < ([3.6] — [8.0]) < 2.9. The YSO population is on average 
redder and fainter than the MS population, with most sources redder than ~ 2.8 in [3.6]- 
[8.0]. Examination of IRAC images of sources from the larger FOV in this YSO CMD region 
suggests that many of them are background galaxies; a dense population of these background 
sources falls within the same 2.5 < ([3.6] — [8.0]) < 2.9 range. 



The CCD (Figure 4(b) ) illustrates that the YSOs are significantly redder than both the 



cluster and the background populations. This is particularly true in [4.5]- [5.8] where only 



- 14 - 



two YSO candidates have colors [4.5] — [5.8] < 0.5. Twenty-three of the 26 YSO candidates 
with these three bands of data he in two sections of color-color space, outlined in green and 
discussed as indicators of strong and weak emission from polycyclic aromatic hydrocarbons 
(PAHs) in Section 14.41 Non-YSOs are strongly clustered around zero on both axes, with no 
non-YSO cluster sources and very few background ("outside cluster") sources falling within 
the YSO color-color spaces. Zero non-YSOs have 1.5 < [4.5] — [5.8]. 



4. Combining Optical and IR: Fitting Spectral Energy Distributions 

We apply the SED fitting tool to a total of 77 "fitter sources," two of which are multiple 
sources with sufficient wavelength coverage to be fit as two separate sources (Sections 14.7.11 
and 14.7.2"]) . Of these 79 sources then, forty-one are YSO candidates, twenty-eight are stellar, 
eight are background galaxies, and two are of unknown nature. Fluxes used are given in 
Table m Images of all YSO candidates are shown in Figure [S] and cover optical, IRAC, 
and MIPS 24 /xm data. 

We characterize 37 sources as well-fit YSOs. We name these with their IRAC pho- 
tometric IDs (used in Table [3]) preceded by Y. The evolutionary stages for the YSOs are 
concentrated in the less evolved Stage I (21 I YSOs) and Stage I or ll (6 l/ll YSOs) with 
a few Stage ll (7 ll YSOs). Of the three remaining well-fit YSO candidates, two may be 
Stage II or Stage ill, and one can be fit with YSO models of all evolutionary stages and is 
considered unclassified. Sources Y090, Y227, Y240, and Y251 are probable YSOs based on 
color, environment, morphology, and best fits. They are Stages I (Y090, Y240, Y251) and ll 
(Y227), but their SED fits are not as good (xVP^^mm > 3-00). 

Table [5] gives the average fitter results for all YSO candidates, including the number 
of optical sources, the optical and IR masses, the YSO's luminosity and accretion rate, 
and the estimated evolutionary stage. Where one or more optical sources correspond to 
a YSO candidate, we attempt to estimate the optical mass in addition to the mass given 
by the fitter. We determine the location of optical sources on the CMD and compare to 



evolutionary tracks (Figure 3(b) ). Only a small percentage of the optical counterparts have 
good photometry and fall within the color-magnitude space covered by evolutionary tracks; 
most are too young, too red, or of too small a mass for these estimates. Intrinsic reddening 
will also cause some optical sources to appear to have somewhat lower masses than they do. 
The optical masses we do calculate (summing over the entire IRAC source) are consistent 
as lower limits with mass estimates from SED fits. 



We fit an additional 38 sources as non-YSO candidates. Nineteen are classified as 



- 15 - 



single, naked stars with Kurucz models (lKurucg|l993l ). and we name them by their IRAC 

photometric ID preceded by K. Nine other sources are clearly stellar (named S ) based on 

imaging but do not have conclusive SED fits. Eight sources are visually identified as galaxies 

(G ). Sources U346 and U703 are the only sources we are unable to classify. They fall 

outside the yl CS" optical FOV, are poorly fit as YSOs, and neither their positions nor their 
morphologies provide clues to their natures. 



4.1. The Fitter 



We use the model SED fitting tool described by lRobitaille et al. I ( 120071 ) , which compares 



fluxe s in specific photometric bands to a grid of 200,000 YSO model SEDs ( jRobitaille et al. 



20061). The model grid c o nsists of model SEDs computed using the radiation transfer codes of 
Whitney et al. I ( l2003al Jbl. l2004j ) for 20,000 sets of physical parameters, 10 viewing angle (from 
edge-on to pole-on), and in 50 different photometric apertures (from 100 to 100,000 AU). 
Among the physical parameters varied in the production of model SEDs are stellar mass, 
radius, and temperature, envelope mass and accretion rate, and disk mass, flaring angle, and 
accretion rate. Ranges for these parameters are determined from observational data. 

A limited number of naked photosphere, asymptotic giant branch (AGB), and galaxy 
SEDs are also fit to the photometric data for comparison . All sources are allowed to be 
fit with YSO model SEDs and Kurucz naked star SEDs ( iKuruczl 119931 ) . including stellar 
photospheres with extinction < Ay < 15 in addition to internal reddening for YSOs. For 
AGBs, we consider extinction in the range 2 < Ay < 20, while galaxy extinction factors are 
3< Ay < 30. 

Synthetic YSO SEDs are computed for five distances in the prescribed range of 60- 
65 kpc, appropriate for the SMC. At a distance of 60.6 kpc, our aperture radii of 1'.'6 in 
IR and optical and 10" in MIPS 24 fim correspond to 9.48 x 10^ AU and 5.92 x 10^ AU 
respectively. The largest aperture available for the computed SEDs is 10^ AU, so this is the 
aperture radius used for our fits. 

The current models are calculated for Solar metallicities and do not include effects from 
PAH emission, external illumination, or protoclusters. Partially as a consequence of the 
physical aperture size at SMC distances, there are often multiple optical sources within the 
aperture of a single IR source. Emission from PAHs can contribute significantly the fiux in 
some IRAC bands. Metallicity is expected to play a role in the timescales of star formation, 
because metals help carry heat away from collapsing gas and dust, (possibly) increasing the 
speed at which a star can form. Models are being developed to include these effects (cf. 
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Sewilo et al. Il2010l ). To account for multiplicity and PAHs, optical fluxes as well as 5.8, 



8.0 /im fluxes must be treated with care, as discussed in Sections |2. 3. H 14. 7. H and 14.7.21 



4.2. Quality of Fits 

We quantify how well a source is fit by a given model SED by considering the value 



of x^/point (as in iRobitaille et al. 1120071 ). In general, the probability that a given model 
reproduces the input data is p ~ e~^^/^ (or ~ — 21np). We consider p to be a Gaussian 
for photometric data points, giving the standard equation (2). In SED fitting, we 

include some fluxes as upper limits only. Our confidence (c = 1 — p) that the source is not 
brighter than this limit is input to the fitter, values for these points are then calculated 
as in equation (3). 

Xdata = 2^ (2) 



=1 
n 



2 _ ^ / -21n(l-c) ifM, >F, 



upper ifMi<Fi 



(3) 



Where Fj is the fiux at wavelength Aj with error bar cTj describing the gaussian probability 
distribution. Mj is the model flux value (which includes extinction and scaling for distance) 
at Aj, and n is the total number of data points. 

Xtotal Xdata ~^ Xupper 

X'^/Pt = xLai/n (i-e-, reduced x^) 

Based on visual examination of both SEDs and images, sources with at least one fit of 
x'^/pt < 3.00 are considered "well-fit" . We average the parameters of all models with x^ /pt < 
i}d IP^min + 0.50) to characterize each YSO candidate. Table [5] lists the number of models 
considered. In all example YSO SEDs (Figures [6], d El |9l [TOl [TTl and [12]), we show the best 
fit (that with Xmin) ^ dark solid line with grey lines indicating all fits considered. SEDs 
of non-YSOs (Figures [TT], [T31 and are shown with only the top twenty-five fits (which 
are representative of the whole set). 



4.3. Characterizing YSOs 



After applying the fitter to all 77 fitter sources, we select YSO candidates. To calculate 
their physical characteristics, we take the well-fit-model parameters and apply a dust-to-gas 
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ratio of 1/6 Galactic, appropriate for star- forming clusters in the SMC (iBot et al.l 120071 ). 



YSO sources can then be classified according to evolutionary stage, based on the frac- 
tiona l disk mass and the env elope accretion rate. We use the physical Stages I, ll, and ill 
from iRobitaille et al. I (120061) wh i ch ar e roughly equivalent to the the observational Class I, 
II, and III classifications of iLadal (119871 ). 



1. Stage l: Embedded Source (Me™ > lO'^M^ yr"^) 

2. Stage ll: Source with Disk (M^™ < lO^^Af^ yr^^ & Mdisk/M^ > 10^^) 

3. Stage ill: Source with Optically Thin Disk or No Disk {M^^v < W'^M^ yr'^ & 
Mrf,,,/M, < 10-6) 



In Table 0, we list the parameters resultant from both the best fit model and the average of 
all models in the x^/pt < ix^/P'^min + 0.50) range. 



4.4. Accounting for Polycyclic Aromatic Hydrocarbons 

Emission from PAHs is not incorporated into the SED models in the current version of 
the fitter. IR flux can be significantly increased b y PAH emission. Of the IRAC bands, only 



[4.5] should be unaffected by PAH contributions (IChurchwell et al. 112004 ). The sources af- 
fected by these PAH emissions are readily identifiable by the characteristic dip in XlogioiRvLx) 
at ~ 4.5/im in the SED, which can be quantified using color-selection such as that shown in 



Figure 4(b) , We must weight the uncontaminated [4.5] band flux more heavily than the mea- 
sured fluxes in the other three IRAC bands. Most appropriately, the heavily contaminated 
[5.8] and [8.0] fluxes would be treated as upper limits (Robitaille, Private communications, 
2009); however, in the absence of longer wavelength measurements, this is impractical. We 
therefore adjust the error bars to 20, 10, 30, and 40% in the 3.6 /xm, 4.5 /im, 5.8 /im, and 
8.0 fim bands, respectively. 

We perform these PAH corrections on the 13 sources falling within the "Strong PAHs" 



area given in Figure 4(b) The relatively blue slope of 3.6 to 4.5 /im, combined with the 
relatively red slope of 4.5 to 5.8 /im, is indicative of a strong dip in the SED at 4.5 fim 
and thus of significant PAH "contamination" in the other IRAC bands. We note sources to 
which we apply PAH corrections in Table O 

Five sources with probable significant PAH contributions (Y217, Y227, Y270, Y287, 
and Y340) do have 24 /im flux measurements; we consider both methods of PAH correction. 
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We perform SED fits considering 5.8 /im and 8.0 fim fluxes as upper limits, relying upon 
the additional SED constraint of the longer wavelength measurements. These are the fits 
we use in our calculation of the mass function and in the determination of YSO stages as 
reported in all figures. The difference in fitting PAHs with 5.8 /im and 8.0 fim fluxes as 
upper limits rather than as having error bars of 30% and 40% is shown in Table |5l For 
sources Y217, Y227, and Y287, fits with increased error bars suggest higher masses and later 
stages of evolution in comparison to fits with 5.8 and 8.0 /im treated as upper hmits. Source 
Y340 is unclassified using the upper limits method but classified as Stage ill using increased 
error bars; the estimated mass is higher using the error bar method. Y270i (Section I4.7.2p 
is the only source which breaks the trend and is fit as Stage I with the error bar method but 
Stage l/ll using the upper limits method. Results using the error bar method are reported 
in Table E] as Y e. 



4.5. YSOs without Optical Counterparts 



Sources Y142, Y143, Y174, Y179, and Y217 (Figure 5(d), 5(e), 5(h), 5(i), and 5(k) 



respectively) have no readily identifiable optical counterparts. Y217, which lies at the tip of 
one of the pillars along the western molecular ridge, is classified as a Stage I or ll source. 
The rest are fit as Stage I and correspond to distinctive dust features. Each is an optically 
dark bump along a molecular ridge. Y142 lies within the "mountain of creation" south of the 
central cluster. All five have masses estimated between 6 and 7 Mq. Y149, Y163, Y493, and 
Y700 lie outside of the ACS FOV and have no known optical counterparts for this reason. 

The SEDs of Y142, Y217, and Y700 are shown in Figures E] and [7] and are markedly 
different. The fit to Y142 is poorly constrained at wavelengths greater than 8 fim and may 
peak anywhere between 8 and 70 /im. Y217 has a double peak SED but is unconstrained 
at wavelengths greater than 24 /xm. Y700's SED is remarkably flat from about 1 /xm to at 
least 24 /xm. 



4.6. YSOs with Single Optical Counterparts 



Sources Y118, Y148, Y251, and Y283 (Figure 5(c), 5(f), 5(m), and 5(o)) correspond to 



single optical point sources, but each is in a different environment. Y118 matches a relatively 
isolated PMS candidate in a region of apparently thin diffuse dust and is fit as a Stage I YSO 
(Figure 8(a)). The mass of the optical source (~ 0.8 Mq) as determined from the optical 
CMD (Figure 3(b) ) is significantly smaller than the ~ 7.4 Mq of the YSO fit, likely indicating 
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the presence of a separate IR-bright source. It is also possible that intrinsic reddening causes 
and underestimate of the optical source's mass; with Ay ~ 3, its mass could be as high as 
1.5 M0, still leaving ~ 6Mq unaccounted for. Stage I source Y148 appears as a very faint 
optical source in a region of optically thick dust. Y251 is poorly fit as a Stage I YSO but 
corresponds to a PMS candidate in a dust peak near cluster center. Y283 is a well-fit Stage ll 
source with an SED peaking around 10 fim (Figure 8(b) ) and an estimated mass of ~ 4.6 M© 
in fair agreement with its optical mass of 3 — 4 Mq . 



4.7. YSOs with Multiple Optical Components 



Approximately 70% of our YSO candidates have multiple optical counterparts. Optical 
sources corresponding to YSO candidates are marked on the CMDs in Figure [31 Most of 
these are fainter or redder than the 1 Myr isochrone, indicating that they are very young 
and probably embedded, as we would predict for sources related to the earliest stages of star 
formation. Where possible we determine the masses of YSOs' optical counterparts through 



comparison the evolutionary tracks (Figure 3(b)); we sum the apparent masses of these 
optical sources to determine the optical masses (Mopt) given in Table [51 We find that YSO 
model masses are systematically higher that optical masses. Many masses are probably 
underestimated, because their embedded nature makes them look fainter in the optical. 
Looking at Figure \5\ there are far more optical sources corresponding to YSO candidates 
than appear on the optical CMD. In particular, sources Y090, Y096, Y223, Y270, Y326, and 



Y327 (Figure 5(a), 5(b), 5(1), 5(m), 5(s), and 5(t)) have more than ten optical counterparts 
each (see also Table [5]). Some of these lie redward of evolutionary tracks on CMDs, preventing 
mass estimates. Others are non-photometric as a result of confusion from source multiplicity 
and high background, as well as their faint nature, and they are not in our catalog (Table [2]). 
There may be embedded IR-dominant sources, likely representing ongoing star-formation 
(as Y340 and Y270), as massive YSOs are expected to evolve more quickly than their low- 
mass counterparts. It is also possible that multiple optical sources are heating more of the 
surrounding ISM than is truly circumstellar matter, causing additional infrared excesses. 
SEDs of distinct multiple sources Y287, Y288, and Y312 are shown in Figures [9l and [TOl 



4.7.1. Source 340 

Optical imaging of Source 340 reveals exactly two optical sources, one optically bright, 
the other faint (Figure [TT]) . They lie at the tip of a dusty pillar highlighted by Ha emission. 
Source 340 is also a PAH source with 24 fim flux. Fitting all data points with either PAH 
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correction method results in a poor Stage III fit (A340; xVP^m ~ 3.5). Fitting different 
wavelength regimes separately results in two good fits, which add together to describe the 
SED nicely. A main sequence fit (K340; X^/ptmin ~ 0-03) accounts for the majority of the 
optical and NIR fiux, while the longer wavelengths are well-fit as an unclassified YSO (Y340; 
X^/ptmin = 0-62) • As with Y217, Y227, and Y287, treating 5.8 and 8 /zm fiuxes as upper 
limits for Y340 results in a lower mass estimate than if we treat them as concrete data points. 
We show the SEDs for A340, K340, Y340, and all three overlaid, demonstrating the improved 
(and more physically meaningful) fit that results from treating K340 and Y340 separately. 
The brighter of the two optical sources distinguishable in Figure 11(b) , has a mass of ~ 4 Mq 
based on its position on the optical CMD (Figure 3(b) ) and likely corresponds to the K340 
fit. The fainter optical source may be dominant in the longer wavelengths, and it is probable 
that one or more protostars are still embedded in the dusty pillar. Physical characterizations 
of each fitting method are given in Table O 



4.7.2. Source 270 

Optically, Source 270 appears to be a protocluster emerging from a pillar of dust (Fig- 
ure [12]). At least ten source s are d iscernible in the optical image, including one B2 star 



classified by iHutchings et al.l (1 19911 . star 12). We fit source 270 in the three ways shown 
in Figure [12 (A270) fitting all data points, (Y270o) fitting optical plus NIR, and (Y270i) 
fitting IRAC plus MIPS 24 /xm. A270 is poorly fit as Stage ll (xVP^^mm = 3-36), while Y270o 
and Y270i are well-fit as Stages I and l/ll YSOs of mass ~ 5.6 Mq (x^/p^mm ~ 1-79) and 
~ 9.3 Mq (x^/P^rrtm ~ 1-05), respectively. Using the error bar method of PAH correction for 
Y270i results in evolutionary Stage I and a slightly lower mass of 8.2 Mq compared to the 
upper limit method. 

Eight of the optical sources corresponding to source 270 have good photometry, and we 



are able to approximate masses for five of these (see Figure 3(b)). Two lie near the MS and 
have masses greater than ~ 5.5 Mq; a third falls on the PMS 3 Mq evolutionary track and 
is much younger /redder than 1 Myr, and the other two are PMS stars with probable masses 
< 1 Mq. Other optical sources are too red or too faint for a mass determination. The sum 
of these optical masses is ~ 15 Mq, which is also the sum of the Y270o and Y270i masses 
given by the fitter. 
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4.8. Non-YSOs 



4.8.1. Stellar Fits 

We are able to fit 18 of the 77 fitter sources (plus K340) as naked stars (MS or G iant 
stars). Four of these correspond to O or B type stars defined by iHutchings et al.l (Il99l[ ) on 
the basis of optical and UV spectroscopy. Five more sources are readily identified as stellar 
and are fit o nly as naked stars but h ave x^/pt ^ 3.00. S293 is one of these and is classified as 
a Bl star by IHutchings et al.l (Il99l[ ). They also report one B2 star corresponding to source 
270 (see Section l4.7.2p . one 06-08 star corresponding to S235, and three B stars that have 
insufficient IR photometry for our SED fitting. A typical stellar SED, for K348, is shown in 



Figure 13(a 



Nine readily identifiable stellar sources (S049, S213, S235, S293, S394, S406, S411, S456, 
and S486) can be decisively fit. S049, S411, and S486 lie outside the ACS FOV. S049 and 
S486 display characteristic diffraction patterns in IRAC imaging and are fit by naked stellar 
templates with 3.0 < X^/ptmin ^ S411 is well-fit but may be fit as a naked star, a YSO, 
or an AGB. The rest fall within the ^C^FOV and are clearly stellar. S293, S394, and S406 
have typical but imperfect photospheric SEDs. 

Three stellar sources are poorly fit because of their evident mid-IR excesses. S456 lies 
on the edge of the ACS FOV in Figure [T6|) and has no reliable optical photometry. Its SED 
does not drop as quickly in the IR as would be expected for a naked star . S213 lies near 



cluste r center, and several optical sources fall within the IRAC aperture. IHutchings et al. 



( 1199 ll ) identify the brightest of these as an 09 star, though there is significant confusion 
even in our high-resolution optical imaging. Its SED falls fairly smoothly through [5.8] but 
then rises sharply th rough 8.0 to 24 /i i n. S2 35 lies at the bright tip of the MS on our optical 



CMD (Figure 3(b) ). IHutchings et al.l ( 119911 ) report that it is (optically) the brightest source 



in NGC 602 and classify it as an 06 or 08 type bright, young MS star based on optical and 
UV spectroscopy respectively. Spectra are shown in their figures 3, 4, and 5 for their star 8. 
Looking at this source in IR wavelengths shows something different. The SED drops steadily 



through [8.0], but the MIPS 24 /im flux is far higher than would be expected (Figure 13(b) ). 



The flux points for optical through IRAC can be well-flt as a MS star, but the 8.0 /im 
measurement is slightly brighter than MS flts would suggest. This is likely the result of the 
O star heating the nearby interstellar medium. Alternate l y, this star could be the flrst to 
emer ge from an embedded cluster (cf. IWalborn et al.lll999l . l2002l : iHeyderi-Malayeri fc Seller 
2OIOI . for similar examples in 30 Doradus and N 66) the O star may exhibit a rn id-IR excess 
due to free-free emission in the stellar wind (as described in iBonanos et al.ll2009l ). or it could 
be a chance superposition. 
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4-8.2. Galaxies and Unidentified Sources 

We visually identify 41 of the 497 IRAC sources with at least one band of data (Table [3]) 
as galaxies in ACS images. Eight of these are fitter sources, six of which are readily identi- 
fiable as ellipticals or spirals and one as an irregular. The eighth galaxy (G211) is outside 
of the ACS FOV but appears markedly elongated in mosaiced IRAC images. All eight are 
best fit as YSOs, some very well. Only seven of the eight are plotted in the IRAC CMD 



(Figure 4(a) ); G150 has ~ 1.5 mag error estimate in [8.0]. Only two of the eight have 5.8 /im 
flux measurements. In Figure [TU we show the SEDs for the two most obvious galaxies in 
the regions, elliptical G133 and face-on grand design spiral G372. These are the only two 
that are not well-fit as YSOs, and we show their best-fit galaxy templates. 



5. THE YSO MASS FUNCTION AND STAR FORMATION RATE 

The total mass of our well-fit YSOs is ~ 300 Mq. We plot a histogram of YSO masses 
(Figure [151 cilso Table [5]) and assume approximate completeness at the peak of the distri- 
bution (~ 6 — 10 Mq) to estimate a lower limit for the total mass. We consider a two-part 

t nr > 0.5 Mq, w here ^ is the 



MF with ^ oc M~t]^j^ for M.^^ < 0.5 Mq a nd £ oc M;,^;' fo r M. 



number of sources with mass Mgtar (as in IWhitney et al.l 12008 



Kroupall200ll . for complete 



IMF discussion). Requiring the two functions to be equal at Mgtar = 0.5 Mq and scaling 
to the peak of the observed mass function, we integrate over the mass range 0.08 — 50 Mq, 
determining a total YSO mass of ~ 2250 Mq. We may then estimate a star formation rate 
(SFR) assuming constant star formation over a given time and considering the physical size 
~ 9.16 X 10"'^ kpc^ of the studied 3' radius region. We consider two possible time scales: 



1. 1 Myr: As seen in Figure [3] (b) the PMS stars related to YSOs are characteristically 
younger than ~ 1 Myr. Applying this time scale, we calculate an SFR of ~ 2.2 x 
10-3 Mq yr-1 (~ 0.24 Mq yr'^ kpc'^). 



2. 2 X 10^ yr: IWhitney et al.l (120081 ) apply a time of 2 x 10^ yr for the formation of current 
YSOs in the LMC, based on the re l ative n umber of Stage I YSOs in their sample and 
the timescale calculations of iLadal (119991 ). For NGC 602, this results in an SFR of 
~ 11.1 X 10-3 Mq yr-1 (~ 1.22 Mq yr'^ kpc'^). 



In our earlier paper (ICignoni et al.l l2009l ). we calculated an SFR for the observed 
optical population of 3 — 7 x 10~^ Mq yr~^. This estimate considers the mass range 
0.45 — 120 Mq and covers only the ACS FOV (the white square in Figure [2]), which has 
a physical size of 58 pc x 58 pc (3.36 x 10^^ kpc^). This results in an optical SFR of 



- 23 - 



0.089 — 0.208 Mq yr~^ kpc~^, assuming star formation has been constant over the past 
2.5 Myr. In order to directly compare this optical estimate with our calculated YSO SFR, 
we integrate the IR MF over the same 0.45 — 120 Mq mass range and obtain a total YSO 



mass of ~ 1300 Mq. ' 
the optical val ue from 



his results in a YSO SFR of 0.14 — 0.71 Mq yr ^ kpc ^. Comparing 



Cignoni et al. 



( 120091 ) (as well as with the similar optical mass function 
determined by lSchmalzl et al.l (j2008l )) with this IR calculation implies that the region's SFR 
has been approximately constant from ~ 2.5 Myr ago to the present, possibly with a slight 
increase in the last 0.5 Myr. 

As a point of comparison, the Orion Nebula is of a similar physical size and age to 
NGC 602. Orion's physical extent is ~ 25 pc radius (cf. iHillenbrandl 119971 ) : in our current 
study, we consider an area of radius ~ 53 pc around NGC 602, more than encompassing 
the main cluster. We esti mate NGC 602 you ng stellar populations to have ages in the range 
~ 10^ yr to a few x 10^ yr. IHillenbrandl (119971 ) quote exactly the same age range for Orion and 
see a similar PMS population. In the same paper, the SFR is calculated as > 10~^ Mq yr~^ 
over the central active star forming area of the Orion Nebula, considering optical sources 
down to a mass of ~ 1 Mq within a FOV of ~ 36' x 34'. If we assume a distance of ~ 400 pc 



Hirota et al.ll2007l ). the physical area is ~ 4 pc x 4 pc or 1.6 x 10 ^ kpc^, and the SFR 



(.e.g. 

is ~ 6.25 MQyr-ikpc-2. Our SFR of 0.14 - 0.71 Mq yr^^ kpc'^ for NGC 602 is a factor of 
ten to fifty less than that given for Orion, unsurprising considering the inclusion of the more 
diffuse star-forming environs of NGC 602's outskirts and the higher molecular gas densities 
in Orion. 



6. SPATIAL AND TEMPORAL DISTRIBUTION OF YOUNG STELLAR 

POPULATIONS 

6.1. Spatial Distribution of Optical and Infrared Populations 

Inspection of the optical CMD (Figure [3]) together with our analysis of combined optical 
and infrared data reveals that different stellar populations coexist in the area: 

1. Young Stellar Objects. YSOs are the least evolved stellar sources. We identify Stage I 
and Stage ll sources with masses 2.5 < M/Mq < 13 (except for Y206 with M/Mq ~ 
26). They are heavily concentrated in the periphery of the cluster, along dusty ridges. 
Most correspond to multiple optical sources with ages less than 1 Myr. 

2. Pre-Main Sequence Stars. Low mass (0.6 < M/Mq < 3) PMS represent the most 
remarkable feature in the optical CMD. They are characterized by red colors (^«f555w~ 
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'^F8i4W ^ 1-0) and faint magnitudes (mF555w ^ 22.3). Their ages are generally less 
than ~5 Myr, certainly less than 10 Myr. They are concentrated near cluster center, 
but clumps also appear in the dusty outskirts of the cluster, corresponding to YSO 
candidates. A total of 494 PMS candidates are included in our final optical photometry 
list (Table ED. 

3. Young Stars. A bright (14.5 < mpsssw ^ 22.3), blue (mpsssw — "^fsmw ^ 0.3) and 
well-defined MS is visible in the upper left of the optical CMD. The majority of these 
stars belong to NGC 602, and the population is strongly concentrated toward the center 
of the cluster. 

4. Old stars. Old SMC wing stars populate the lower part of the MS (mpsssw ^ 22.0). A 
few red giant branch (RGB) and red clump stars are visible in the upper right of the 
CMD. These are primarily a background population. They are m ost visible around the 



edges of our analyzed region, outside the nebula. Our work in ICignoni et al.l (120091 ) 
indicates that this population is consistent with the nearby field population and at a 
larger distance modulus than the young population. 

Figure [TB] outlines the spatial distribution of these stellar populations with a grid on 
the ACS Ha. PMS stars are marked as blue dots and YSOs as circles with colors indicating 
their determined stages. Stellar sources are marked as squares. Background galaxies are 
marked as diamonds. Figure [17] shows optical CMDs corresponding to each of the grid 



regions delineated in Figure [16] with isochrones from Figure 3(a) for reference. The number 



of Spitzer-identiiied YSO candidates is noted in the top right-hand corner of each CMD. 



PMS candidates associated with Spitzer YSOs are shown in red (as in Figure 3(a)). 



The concentration of PMS stars versus YSOs is striking. Of the candidate PMS stars, 
38% lie within grid-section (4,3). The majority of YSO candidates lie farther from cluster 
center, notably along ridges in the East/South-East and West, and four YSOs (Y149, Y163, 
Y493, and Y700) lie outside the optical FOV. Grid section (5,3) contains 6 of the 36 well- 
fit YSOs but less than 7% of the PMS stars. Grid-section (4,2) includes 6 well-fit YSOs 
(seven if Y270o and Y270i are counted separately), 3 YSO candidates that are not well- 
fit, and another probable YSO with insufficient photometry for our fitting requirements. 
Although these (4,2) YSOs correspond to many PMS candidates, less than 12% of all PMS 
candidates fall within this grid-section. Also notable are the evident ages of the PMS sources 
in the central versus outer grid- sect ions. In (4,3), the PMS population aligns well with the 
1 Myr isochrone (red) with some older and some younger. Other grid-sections (e.g., (3,2), 
(3,3), (3,4), (4,2), (4,4), (4,5), (5,3), and (6,3)) suggest younger PMS populations that lie 
redward of the 1 Myr isochrone (in red). While some of this is likely the result of increased 
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reddening in these dustier areas, extinction is insufficient to account for the measurable 
trend. (Figure 3(a) shows the appropriate reddening vector.) 

We consider two components of the MS, old and young populations. Many of the edge 
regions contain only old MS a nd RGBs. Vir t ually all of the optical sources in grid-sections 
(1,3) and (1,5), which contain Icignoni et all fcoooh sub-clust ers NGC 602-B and NGC 602- 



Schmalzl et al.ll2008l ). are old stars. The young MS 



B2 (alternately B 164 and Cluster A in 
(mi7555vi/ < 17) population is concentrated at cluster center, particularly in grid-section (4,3) 
with the PMS concentration. 



6.2. Progression of Star Formation in NGC 602 

The distances from cluster center to PMS and YSO candidates of different stages reveal 
a gradient. We define the cluster center as the center of the distribution of PMS candi- 
dates (RA, Dec: 22°. 3855, -73°.5584 or l'^29™32^5, -73°33'30'.'24), which is approximately 
the center of the main 0/B association. We measure the projected distance to PMS and 
YSO candidates. As a correction for the three-dimensional distance, we make the simple 
assumption that the actual physical distances are a/2 times these projected distances. 

We have a list of 565 PMS candidates, 484 of which lie within 30 pc of cluster center. 
We consider these 484 to be strong candidates, as their positions are consistent with a 
common triggering process, and their clustering near the 0/B association makes their nature 
more certain. The primary locus of the PMS population is clear from image examination 
(See Figures and [TS]), and the distribution is approximately gaussian. We define the 
PMS distance as the radius from cluster center that will include 50% of these strong PMS 
candidates and use error bars defined by the radii including 40% and 60% of the strong 
candidates. The result is Tppjs = 9.0 ± f l pc. 

For the average distances to YSOs, we remove sources outside the ACS FOV because 
their nature is less certain, and their formation may be triggered in a different manner. Y700, 
for example, lies on its own semi-circular ridge of 8.0 fim emission. The average distances 
from cluster center to YSO stages are: Stage ll ~ 12.0 ± 5.5 pc. Stage l/ll ~ 16.6 ± 5.9 pc, 
and Stage I ~ 19.4 ± 9.2 pc, where the error bars are the standard deviations of the mean. 

We consider a simple picture of NGC 602 in which the region consists of a central O /B 
association surrounded by concentric shells of star formation. We make the assumption that 
star formation began at cluster center (the center of the PMS distribution) approximately 
3 Myr ago and is propagating outward. This picture is shown schematically in Figure [181 for 
simplicity, we only include shells for PMS, Stage ll, and Stage I sources. As mentioned in 
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Section [Tj lNigra et al.l (|2008[ ) have shown that the nebular expansion velocities are negligible, 
so collapse/star formation (beginning YSO Stage 0) is most likely being triggered along 
ionization fronts. We assume a local sound speed of ~ 10 km/s (as appropriate for an 
ionized region at T~ lO^'K) for the speed of star formation propagation. Dividing the 
distances by this sound speed, we derive timescales for the populations. The PMS population 
formed within ~ 0.9 Myr of the 0/B association formation. The average time between the 
formation of the 0/B association and the YSOs increase with less evolved YSO stages as 
follows: Stage II ~ 1.2±0.55 Myr, Stage l/ll ~ 1.66 ±5.9 Myr; and Stage I ~ 1.94±0.92 Myr. 
The Stage l/ll shell forms between Stage I and Stage ll, representing a mix of Stage I and ll 
sources. These timescales associated with the YSO stages represent samples of the continuous 
propagation of star formation radially outwards in the NGC 602 region; there is, however, 
overlap between all of the shells. Sources of uncertainty include imprecise distances due to 
three-dimensional structure and low-number statistics. 

The relative times between the triggering of these Stages can also be used to estimate 
Stage lifetimes. PMS stars and what are now Stage ll sources are triggered 3 ± 6 x 10^ yr 
apart, indicating that the Stage II evolution lasts on the order of ~ 3 x 10^ yr. Likewise, 
Stage I sources have a lifetime ~ 3 x 10^ yr, from the time difference between Stage I and ll 
shells. We are unable to account for the probable effects of mass on the range of YSO evo- 
lutionary timescales. The smallest value for the Stage I shell radius (i.e., the best estimate 
minus the uncertainty) is ~ 10 pc (1 Myr), indicating that triggering occurred as much as 
2 Myr, which may be taken as a firm upper limit to Stage I lifetimes. Massive sources are 
expected to evolve more quickly than low mass sources; most of our candidates are protoclus- 
ters/multiple sources, leading us to expect a multiplicity of masses and thus stages within 
individual YSO candidate. The uncertainty in these numbers is large, particularly because 
we are unable to analyze the expected correlation between mass and stage lifetime with such 
low number statistics and the current resolution. However, the results are interesting; the 
temporal differences between the triggering of PMS and Stage ll. Stage ll and Stage l/ll, 
and Stage l/ll and Stage I are all approximately equal (a few 10^ yr). Our Stage I lifetime 
is consistent with the li fetim es measured for YSOs in the galactic star-forming region, M17, 
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We construct a time-ordered scenario for the progression of star formation in NGC 602. 
Approximately 7 Myr ago, two expanding Hi shells began to interact, creating an over- 
density in the region that is now N90. After ~ 3 Myr, turbulenc e subsided sufficiently for 



star formation to begin in earnest, ~ 4 Myr ago (INigra et al.l 120081 ). The current population 
of bright main sequence (O and B) stars and low-mass PMS stars formed near the center 
of the over-dense region approximately 2-3 Myr ago. The radiation from the massive stars 
began to erode the surrounding nebula, creating a photodissociation region and triggering 
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further star formation around the edges. The formation of the current Stage ll YSOs was 
triggered ~ 1.8 Myr ago, the formation of the current Stage I YSOs ~ 0.7 Myr later (or 
~ 1.1 Myr). 



7. CONCLUSIONS 

We present a muhi-wavelength analysis of photometry and imaging of the NGC 602 
active star-forming region in the SMC, covering HST 0.55 through Spitzer 24 fim. From 
these data, we define stellar and proto-stellar populations and their spatial distribution. We 
estimate the present-day star formation rate and derive time scales for the formation of 
Stage I and Stage ll YSOs. We provide full mutli-wavelength photometric catalogs online 
and present approximate physical and evolutionary parameters for all of our YSO candidates. 

Our primary modes of source analysis are CMD examination and SED fitting. Optical 
HST CMDs reveal ~565 PMS candidates, essentially low mass Stage ill YSOs, through 
isochrone fitting. Through multi-wavelength SED fitting, we identify 41 YSO candidates, 
including 24 Stage I, 8 Stage l/ll, 5 Stage ll, 2 Stage ii/iii, and 2 unclassified candidates. 
High-resoultion HST imaging shows that > 70% of the YSO candidates include multiple 
sources or are protoclusters, and most of these optical sources are PMS candidates. Efforts 
to construct YSO protocluster models and incorporate them into the SED fitter are underway 
but are beyond the scope of this paper. For the ~20 YSO candidates, we are able to estimate 
masses for one or more optical counter-parts via comparison with CMD evolutionary tracks, 
and we find consistency between these lower limit optical masses and the YSO SED fitter 
masses. We also construct a mass function from YSO SED fitter masses and derive a present- 
day star formation rate of 0.2-1 Mq yr~^ kpc~^. 

Finally, we present a quantitative analysis of the spatial distribution of the YSO popu- 
lation with respect to the central cluster and PMS population. We find that star formation 
has progressed from cluster center to the edge of the star forming dust cloud in NGC 602. 
The PMS stars are heavily concentrated near cluster center and that the YSO population 
distribution can be represented as concentric shells with Stage ll sources p referentially closer 



to cluster center and Stage I sources farther away. Previous observations (INigra et al.ll2008[ ) 
have shown that there is no significant expansion of the dust shell in which most of our YSO 
candidates are located and that the photo-ionization front is the prime mover of the star 
formation activity. We therefore correlate average distances of the Stage I and ll YSOs from 
cluster center with the times at which their formation is triggered; we divide the distances 
by the sound speed. Relating the timescales, we find the lifetimes of each YSO Stage to be 
a few 10^ yr, comparable to timescale estimates in the literature which apply independent 
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techniques to galactic star-formation regions. 
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Fig. 1. — Image of NGC 602 in optical through 24 /im. Circle of radii 1'.'6 indicate the 
positions of YSO candidates, color-coded by evolutionary stage as in the legend. The addition 
of a white circle indicates that we do not consider the fit good. North is up, east to the left. 



- 33 - 




RA ■ . 1. 

(J2000) 22.8 22.6 22.4 22.2 . 22.0 • 21.81 

_■ i—l , , , I , i—i I , , , I , , , J , , !_■ 



Fig. 2. — Large 8.0 fim image of the NGC 602 region. Black Box: Deep Spitzer/IRAC FOV; 
White Box: HST/ACS/WFCFOV; Large White Circle: 3' radius region studied here. Small 
white markings: 77 "fitter sources." The full area shown is 18' x 18'. Sources outside the 
Cluster Region are marked in grey in Figure H] to represent the background population. 
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Fig. 3. — (a) Optical Color-Magnitude Diagram of the entire ACS field of view, using 
apparent magnitudes in the HST/ACS/WFC hands F555W and F814W (roughly Fand /, 
respectively). The PMS lies to the right (red) of the lower MS and is a few times 10^ yr. 
We also mark optical sources associated with our stellar and proto-stellar fitter sources. 
Blue asterisks are spatially related to IR sources fit as (non-YSO) stars. Red diamonds 
are spatially related to fit YSO sources; many "YSOs" e ncompass mult iple optical sources, 



mostly P MS. The 1 and 10 Myr FRANEC isochrones from lCienoni et all ( Z = 0.004; 120091 ) 



based on lSiess et al.l (120001 ) , are s hown as representative age boundaries. The 6 Gyr isochrone 



is based on lBertelli et al.l (jl994| ) with Z = 0.001. The yellow line shows Equation [H points 
to the right are PMS candidates. The arrow indicates the reddening vector for Ay = 0.24 
[E[B — V) = 0.08). (b) CMD used to determine masses of YSOs' optical counterparts. Blues 
show the evolutionary tracks for various stellar masses from 0.5 to 5.5 Mq. The dashed black 
line shows the 1 Myr isochrone. Colored points are related to fit YSO candidates. 
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Fig. 4. — (a) IRAC CMD with object types. Colored points are for the 3' cluster region with 
YSO candidates marked as dots, color-coded by evolutionary stage. Black dots are sources of 
unknown nature within the 3' region. Grey points are for sources in an area of 18' x 18' outside 
the main 3' radius cluster region, (b) Color-Color diagram with object types and PAH color- 
color space indicated. Boxes outline empirically determined regions with significant PAH 
contribution by selecting sources that are both relatively blue in [3.6]- [4.5] and relatively red 
in [4.5]- [5.8]. The most prominent 4.5 fim dip in the SEDs correspond to the color-color 
region ([3.6] - [4.5]) < 0.5 and 1.5 < ([4.5] - [5.8]). Sources with shallow 4.5 /xm SED 
dips (slight PAH contribution) are found in the color region 0.3 < ([3.6] — [4.5]) < 0.8 and 
0.8 < ([4.5] — [5.8]) < 1.5. (For clarity, we have allowed a slight space between the two boxes 
on the plot.) 
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Fig. 5. — YSOs in NGC 602. Each square image is approximately 15" (4.5 pc) on a 
side. Circles are 1'.'6 in radius, representing the IRAC photometric apertures. Circle col- 
ors are Red (Stage l), Orange (Stage i/n), Yellow (Stage ll), Green (Stage ii/iii, Lavender 
(Unclassified), and White (Poorly fit). The underlying images are 8-color: Blue= Optical 
(F555W[1^ + F658N[Ha] + F814W[/] from HST/ACS). Green to Yellow= Infrared IRAC 
3.6 + 4.5 + 5.8 fim. Red= IRAC 8.0 /im. Magenta - MIPS 24/im. The regions shown in 



5(g) , 5(x) , and 5(y) are outside the ACS FOV. 
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(a) 



(b) 



Fig. 6. — Example SEDs for YSO candidates with no optical counterpart. Both SEDs are fit 



with optical data used as upper limits. Images of the two sources are shown in Figures 5(d) 
and 5(k), (a) SED for Stage I source Y142. The fit indicates a mass of ~ 7.4 Mq. With no 
24 /im measurement, the fit is poorly constrained, (b) SED of Stage l/ll source Y217 with 
mass ~ 5.7 M0. The deep dip at 4.5 yum indicates significant PAH contamination, and in 
the presence of the 24 /im measurements, 5.8 and 8.0 /im fluxes are treated as upper limits. 
The dashed grey line (Y217e) shows the best fit using the error bar method of correcting for 
PAHs and shows quite a different shape from the fit using the upper limit method. Filled 
circles - photometric data points with applied error bars. Triangles - photometric data points 
used as upper limits in SED fitting as described in the text. All fits within 0.5 x^/pt ^'^'^ 
shown in grey with the best fit in black. The dashed black lines indicate the photosphere 
corresponding to the best fit. All SEDs are available online. 
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Fig. 7. — Y700 lies outside the ACS FOV along a ridge of dust east of the main clusters, 
highlighted in 8.0 /zm in Figure [H It is fit as a Stage I source with a flat SED. Based on this 



fit, the source's mass is ~ 4.3 Mq. See also Figure 5(y) Symbols as in Figure [61 
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(a) 

Fig. 8. — Example SEDs for YSO candidates with single optical counterparts. Images are 



Figures 5(c) and 5(o) (a) Y118 a relatively isolated Stage I source with a mass estimate of 
~ 7.4 M0 based on this SED but approximately corresponding to a single PMS star (cyan 



in Figure 3(b) ) with mass ~ 0.8 M0. (b) Y283 is a bright Stage ll source of mass ~ 4.5 Mq; 



it is plotted in yellow in Figure 3(b) Symbols as in Figure El 
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Fig. 9. — Example SEDs for YSO candidates with multiple optical counterparts. Y287 and 
Y288 are fit as Stage ll and I, respectively, and lie near each other on the plane of the sky. 
In Figure 5(q), they appear to lie along adjacent molecular ridges. Corresponding optical 
sources are marked in green in Figure 3(b) Symbols as in Figure El 




(a) 



(b) 



Fig. 10. — Y312 appears as an arc of Pre-Main Sequence stars (also Figure 5(r) , five orange 



optical sources in CMD Figure 3(b) ) fit as a Stage I source. SED symbols as in Figure ID 
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Fig. 11. — Source 340. An optically bright source dominates (a) the three-color optical image, 
but (b) the narrow band Ha reveals a second, fainter optical point source, (c) In eight colors, 
the infrared source looks orange (strong emission in IRAC bands); the approximate photometric 
aperture is shown. Images (a), (b), and (c) are approximately 22" on each side. The three black 
and white SEDs are (d) the K340 stellar optical fit, (e) the Y340 IR fit, and (f) the A340 fit to aU 
photometric data, (g) The color SED shows the three fits overlaid together with A340 in black, 
K340 in blue, and Y340 in violet. SED symbols as in Figure O 



(a) Optical ACS: F555W + 
F658N + F814W 



(b) ACS: Ha 



(c) Eight-color 




Fig. 12. — Images and SEDs for source 270. Images (a), (b), and (c) are appro ximately 15" 
on each side. The brightest optical source is a B2 star ( iHutchings et al.lll99ll ). We fit the 
photometric data best as two Stage I YSOs. The photometric optical sources are plotted in 
red in Figure 3(b) SED symbols as in Figure El 
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Fig. 13. — Two stellar SEDs. K348 is a well-fit naked star. According to iHutchings et al. 



(I1991I ) S235 is an 06 star, optically the brightest in NGC 602. Our data is well fit as a naked 
star if we ignore the MIPS 24 fim flux, which is likely the result of the O star heating the 
nearby interstellar medium. Symbols as in Figure |6l 
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0.1 1 10 100 0.1 1 10 100 

X (/xm) A {/j,m) 



(a) (b) 




(c) (d) 

Fig. 14. — Example SEDs for known galaxies. G133 is a readily identifiable as a bright 
elliptical and appears bright green in the west of Figure [TJ It is fit best (though poorly 
fit) by AGB templates; here, we show the best fit galaxy template. G372 is the face-on 
grand design spiral just north-east of the main nebula in Figure [1] It is best (though not 
well) fit by YSO models; we show the best galaxy template. The other 6 galaxies in our 77 
fitter sources are all well-fit as YSOs. None of these galaxies are well-fit by galaxy templates 
currently incorporated in the fitter. All galaxy SEDs in online figures show the best-fit 
galaxy templates. Symbols as in Figure Ei 
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Fig. 15. — Mass Function of YSO candidates. The total mass of the 41 fit YSO cand idates is 
30 Mq. The total integrated mass, assuming a two-part mass function (grey line; iKroupa 
200lh from 0.08 to 50 Mq, is ~ 2250 Mq. 
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Fig. 16. — Distribution of young sources on an inverted F658N (Ha) ACS image. Circles 
indicate IRAC YSO candidates and are shown with the IRAC photometric aperture of 1'.'6 
radius and are color-coded by object type as in legend. White circles indicate probable YSOs 
(based on their IR colors and environments) that do not have enough bands for SED fits. 
Sources for which we can estimate stages but which are not well-fit are indicated with white 
in ad dition to their appropri ate stage colors. Boxes indicate fit stellar sources and 0/B stars 
from iHutchings et al.l (Il99ll ). All 565 PMS candidates identified in the optical are plotted. 
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012 012 012 012 012 

m(F555W)-m(F814W) 



Fig. 17. — Optical CMDs for grid sections in Figure [161 The youngest optical sources (both 
PMS and upper MS) are concentrated near cluster center. YSOs are concentrated farther 
from cluster center, especially along the "broken ring" of dusty ridges. The older population 
is more diffuse and appears to be largely background. The number of Spitzer-identified YSO 
candidates within each grid region is noted. Optical sources related to YSOs are marked 
in red; optical sources corresponding to fit IR stars are marked in blue. Isochrones are for 
IMyr (dark red), 10 Myr (light green in sections (3,3), (3,4), (4,1), and (4,2)), 50 Myr (true 
green), and 6 Gyr (blue). 
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Fig. 18. — Relative timescales for triggered formation of YSO Stages. We assume a sound 
speed Cs ~ 10 km s~^ and calculate the time from cluster center (the center of the PMS 
distribution) to the average positions of YSO Stages, multiplying by a factor \/2 as a partial 
correction from line-of-sight distance from cluster center. The large cyan circle encloses 50% 
of the PMS distribution within ~ 30 pc of cluster center (blue dots). The red and yellow 
circles show the average radial distances to Stage I and Stage ll YSOs. Background image: 
ACSRa. 
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Table 1. NGC 602 log-book of optical observations. 



Image Name 


Filter 


Exp. Time 


R.A. 


Dec. 


Instrument 






(sec.) 


(J2000) 


(J2000) 




J92F05LIQ 


F555W 


3 


01:29:27.57 


-73:33:17.1 


HST/ACS/WFC 


J92FA6R7Q 


F555W 


3 


01:29:27.57 


-73:33:17.1 


HST/ACS/WFC 


J92F05LJQ 


F555W 


430 


01:29:27.57 


-73:33:17.1 


HST/ACS/WFC 


J92F05MOQ 


F555W 


430 


01:29:27.25 


-73:33:17.4 


HST/ACS/WFC 


J92F05LQQ 


F555W 


430 


01:29:27.89 


-73:33:16.9 


HST/ACS/WFC 


J92F05LSQ 


F555W 


430 


01:29:28.21 


-73:33:16.7 


HST/ACS/WFC 


J92F05LUQ 


F555W 


430 


01:29:28.52 


-73:33:16.4 


HST/ACS/WFC 


J92F05MSQ 


F814W 


2 


01:29:27.57 


-73:33:17.1 


HST/ACS/WFC 


J92FA6QZQ 


F814W 


2 


01:29:27.57 


-73:33:17.1 


HST/ACS/WFC 


J92F05MTQ 


F814W 


453 


01:29:27.25 


-73:33:17.4 


HST/ACS/WFC 


J92F05N9Q 


F814W 


453 


01:29:27.57 


-73:33:17.1 


HST/ACS/WFC 


J92F05N7Q 


F814W 


453 


01:29:27.89 


-73:33:16.9 


HST/ACS/WFC 


J92F05N3Q 


F814W 


453 


01:29:28.21 


-73:33:16.7 


HST/ACS/WFC 


J92F05MXQ 


F814W 


453 


01:29:28.52 


-73:33:16.4 


HST/ACS/WFC 


J92FA6R0Q 


F658N 


636 


01:29:19.35 


-73:33:18.7 


HST/ACS/WFC 


J92FA6R2Q 


F658N 


636 


01:29:18.71 


-73:33:19.0 


HST/ACS/WFC 


J92FA6R4Q 


F658N 


636 


01:29:18.08 


-73:33:19.3 


HST/ACS/WFC 


rccd061014.059 


V 


45 


01:29:24.78 


-73:33:15.6 


SMARTS/ANDICAM 


rccd061014.060 


V 


45 


01:29:24.80 


-73:33:15.6 


SMARTS /ANDICAM 


rccd061014.061 


V 


45 


01:29:24.81 


-73:33:15.4 


SMARTS/ANDICAM 


rccd061014.062 


I 


30 


01:29:24.82 


-73:33:15.5 


SMARTS/ANDICAM 


rccd061014.063 


I 


30 


01:29:24.83 


-73:33:15.7 


SMARTS/ANDICAM 


rccd061014.064 


I 


30 


01:29:24.85 


-73:33:15.5 


SMARTS/ANDICAM 
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Table 2. Optical photometry table. The full 4535 -source table is availa ble online and 
includes matches to IRAC sources, and optical ISchmalzl et al.l ( l2008l ) stars and 



Hutchings et al.l (Il991f l 0/B stars (matched within 0'.'2). The ACS filters F555W and 



F814W are comparable to Fand / , respectively. Flags indicate the photometric data used: 
2 ACS, 6 ground-based (SMARTS), 12 ACS/SMARTS combination. Designations begin 

"HSTNGC602 ." 



ID Designation R.A. Dec. F555W dF555W F814W dF814W Flag 

(°, J2000) (°, J2000) (mag) (mag) (mag) (mag) 



1 


J012935. 538-733330.75 


22.39806 


-73.55852 


17.257 


0.024 


17.310 


0.076 


6 


2 


J012950.969- 733139.21 


22.46245 


-73.52756 






17.060 


0.093 


6 


3 


J012934.757-733247.31 


22.39484 


-73.54645 


17.603 


0.014 


17.594 


0.036 


6 


4 


J012944.118- 733206.11 


22.43388 


-73.53502 


17.695 


0.014 


17.853 


0.089 


6 


5 


J012937.031-733325.52 


22.40429 


-73.55707 


17.742 


0.015 


16.734 


0.014 


12 


6 


J012935. 517-733329.37 


22.39798 


-73.55814 






17.978 


0.029 


6 


7 


J012921. 403- 733251. 69 


22.33920 


-73.54766 


17.897 


0.016 


17.342 


0.095 


6 


8 


J012926. 544- 733230.95 


22.36064 


-73.54190 


18.012 


0.017 


17.714 


0.030 


6 


9 


J012930.638- 733406.40 


22.37761 


-73.56842 


18.069 


0.021 


18.042 


0.026 


6 


20 


J012935. 361-733300.21 


22.39735 


-73.55004 


18.202 


0.038 


17.958 


0.031 


2 
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Table 3. IRAC photometry table. The full 497-source table is available onlin e and notes 
matches to optical sources from Table [2] as well as matches to iGouliermis et al.l (120071 ) YSO 
candidates. Designations begin with the prefix "SSTNGC602 ." 



ID 


Designation 


R.A. 


Dec 


[3.6] 


d[3.6] 


[4.5] 


d[4.5] 


[5.8] 


d[5.8] 


[8.0] 


d[8.0] 






(°, J2000) 


{°, J2000) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


(mag) 


194 


J012920.73- 733327.1 


22.3364 


-73.5575 


15.040 


0.016 


15.027 


0.027 


15.101 


0.108 


15.339 


0.262 


195 


J012938.80- 733411.4 


22.4117 


-73.5698 


16.142 


0.087 


15.926 


0.084 










196 


J012919.88-733322.5 


22.3328 


-73.5563 


16.116 


0.033 


15.768 


0.046 


14.058 


0.044 


12.466 


0.021 


197 


J012918. 95-733319. 4 


22.3290 


-73.5554 


15.927 


0.027 


15.604 


0.037 






12.429 


0.020 


198 


J012936.38- 733403.6 


22.4016 


-73.5677 


14.507 


0.021 


14.182 


0.030 


13.698 


0.047 


12.104 


0.034 


199 


J012933.75- 733357.0 


22.3907 


-73.5658 


17.298 


0.234 


16.053 


0.124 










201 


J012921. 57-733325.4 


22.3400 


-73.5571 


15.040 


0.016 


15.518 


0.039 










203 


J012925.36- 733331.1 


22.3557 


-73.5587 


17.720 


0.116 


17.106 


0.138 










204 


J012922. 91-733324.7 


22.3455 


-73.5569 


17.733 


0.108 


17.247 


0.152 










205 


J012927.67-733334.9 


22.3653 


-73.5597 


15.240 


0.018 


15.264 


0.034 


15.712 


0.174 






206 


J012920.48-733316.8 


22.3354 


-73.5547 


16.097 


0.031 


15.646 


0.039 


13.302 


0.025 


12.584 


0.022 



Table 4. Fluxes of all 77 fitter sources, including those fit in multiple ways. Uncertainties are included in the full 

online photometry table and are on the order of 10 to 20%. 



Name R.A. Dec. Fluxy Fluxj- Fluxj Fluxjj FIuxr-^ Flux3.e^„ Flux4.5j,„, Fluxs.g^m Fluxg^^ Flux24^m 

(°, J2000) (°, J2000) (mjy) (mjy) (mjy) (mjy) (mj) (mjy) (mjy) (mjy) (mjy) (mjy) 



Y090 


22.36823 


-73.67969 


0.0003^^ 


0.0012^^ 


0.017 






0.12 


0.054 


0.43 


1.1 




Y096 


22.27672 


-73.56351 


0.0075^^ 


0.012^'' 


0.189° 


0.35™ 


0.30™ 


0.93 


1.0 


1.9 


4.0 


36 


Y118 


22.37345 


-73.57505 


0.0018 


0.0043 








0.019 


0.044 




0.070 




Y142 


22.36902 


-73.56999 


0.0003^^ 


o.ooos''^ 








0.097 


0.089 


0.14 


0.41 




Y143 


22.31988 


-73.66133 


0.0039 


0.0048 








0.049 


0.045 




0.31 




Y148 


22.39630 


-73.67331 




0.0001 








0.11 


0.075 


0.27 


0.69 




Y149 


22.23401 


-73.64612 






0.036 


0.058 


0.068 


0.066 


0.057 




0.17 




Y162 


22.38200 


-73.66926 


0.0010 


0.0019 


0.032^° 


0.047^0 


0.068™ 


0.33 


0.39 


0.66 


1.6 


6.3 


Y163 


22.24721 


-73.54578 






0.034 


0.047 


0.063 


0.046 


0.042 




0.033 




Y170 


22.37877 


-73.56824 




0.017^^ 


0.022 


0.044 




0.16 


0.15 


0.27 


0.82 




Y171 


22.38497 


-73.56902 


0.0012 


0.0015 


0.026 






0.34 


0.40 


0.66 


1.5 




Y174 


22.39191 


-73.56994 


0.0003'''°' 


0.0008"^ 


1.0™ 


0.99™ 


0.63™ 


0.14 


0.14 


0.21 


0.57 




Y179 


22.32686 


-73.55694 


0.0037 


0.0033 








0.064 


0.052 




0.31 




Y196 


22.33283 


-73.66628 


0.0018 


0.0022 


0.029 


0.033 




0.10 


0.089 


0.27 


0.66 




Y197 


22.32897 


-73.66640 


0.016 


0.0044 








0.12 


0.10 




0.68 




Y198 


22.40160 


-73.56769 


0.037 


0.20 


0.62 


0.69 


0.61 


0.44 


0.38 


0.38 


0.92 




Y2()6 


22.33535 


-73.55467 


0.0028'''°' 


0.0025^^ 


0.028™ 


0.043^° 




0.10 


0.099 


0.55 


0.59 




Y217 


22.33283 


-73.55347 


0.000005'^-'' 


0.00008 








0.19 


0.16 


0.5489 


1.3 


1.5 


Y223 


22.39958 


-73.56438 


0.0019 


0.0053 


0.066™ 


0.090^° 


0.094^0 


0.37 




1.0 


2.5 




Y227 


22.40573 


-73.56458 


0.0043^^ 


0.0039^^ 








0.67 


0.48 


1.7 


4.3 


11 


Y237 


22.33603 


-73.66176 


0.0017^^ 


0.0018^^ 


0.36™ 


0.21™ 


0.14™ 


0.19 


0.20 


0.38 


0.80 


1.9 


Y240 


22.40832 


-73.66469 


0.031 


0.018 








0.67 


0.48 




4.3 




Y251 


22.39863 


-73.66098 


0.0009 


0.0015 


0.012™ 


0.020™ 




0.18 


0.12 


0.67 


1.4 




Y255 


22.40668 


-73.66194 


0.079 


0.066 








0.071 


0.089 




0.32 




Y264 


22.40414 


-73.56005 


0.13 


0.10 








0.16 


0.15 




0.46 




Y270o 


22.39799 


-73.55843 


Q 4795 


().320-' 


0.43''" 


O.45OO 


0.39™ 


0.70 


0.50 


1.7 


4.2 


37 


Y270i 


22.39799 


-73.56843 


0.47 


0.32 


0.43 


0.45 


0.40 


0.70^^ 


0.50^^ 


1.7™ 


4.2™ 


3795 


A270 


22.39799 


-73.66843 


0.47 


0.32 


0.43 


0.45 


0.39 


0.70 


0.50 


1.7 


4.2 


37 


Y271 


22.38932 


-73.66647 


0.016 


0.019 








0.031 


0.030 




0.066 




Y283 


22.37687 


-73.66292 


0.0083 


0.011 


0.13™ 


0.15™ 


0.21™ 


0.17 


0.20 


0.32 


0.67 


1.3 


Y286 


22.34611 


-73.64764 


0.0024 


0.0041 


0.062 


0.082 


0.081 


0.12 


0.097 


0.30 


0.66 




Y287 


22.43070 


-73.56199 


0.0009^** 


0.0013^^ 


0.016 






0.13 


0.12 


0.35 


0.78 


0.70 


Y288 


22.17680 


-73.01200 


0.0005^"^ 


O.OOIO^'^ 




0.023 




0.13 


0.15 


0.34 


0.73 




Y290 


22.40450 


-73.55705 


0.026 


0.042 


0.6830 


0.7490 


0.47™ 


0.28 


0.22 


0.17 


0.31 




Y312 


22.41324 


-73.55526 


0.0006'''°' 


0.0016^'^ 


0.024 


0.033 


0.041 


0.11 


0.077 


0.30 


0.79 




Y326 


22.43648 


-73.55702 


0.0016^^ 


0.0037^^ 


0.023^0 


0.039^0 




0.11 


0.069 


0.22 


0.58 




Y327 


22.39873 


-73.66069 


0.016"^ 


0.019^^ 


0.18™ 


0.18™ 


0.37™ 


1.9 


2.5 


4.6 


8.6 


39 


Y340 


22.36068 


-73.64176 












0.10 


0.060 


0.20 


0.61 


0.66 


A340 


22.36068 


-73.64176 


0.24 


0.22 


0.16 


0.11 


0.082 


0.10 


0.060 


0.20 


0.61 


0.66 


Y368 


22.39666 


-73.64616 


0.0036 


0.0050 








0.016 


0.021 




0.060 




Y396 


22.35085 


-73.53134 


0.042 


0.04 


0.016 






0.060 


0.053 




0.037 




Y493 


22.31646 


-73.50489 






0.036 


0.080 


0.11 


0.089 


0.075 




0.095 




Y700 


22.62856 


-73.54969 






0.020 






0.069 


0.095 




0.24 


0.73 


S049 


22.26367 


-73.67036 






4.0 


41 


28 


13 


8.0 


4.9 


2.9 


0.71 


K050 


22.38160 


-73.69061 






2.2 


2.1 


1.6 


0.71 


0.42 


0.28 


0.13 




K063 


22.27820 


-73.67049 


0.22 


0.41 


0.86 


1.3 


0.83 


0.43 


0.23 


0.14 


0.12 




K086 


22.30000 


-73.67027 


0.20 


0.33 


0.62 


0.64 


0.39 


0.19 


0.12 


0.062 


0.060 




K181 


22.40200 


-73.67003 


1.1 


2.7 


3.9 


5.3 


3.3 


1.7 


1.1 


0.71 


0.74 





Table 4 — Continued 



Name 


R.A. 


Dec. 


Fluxv 


Flux/ 


Fluxj 


Flux/r 


Flux/f 




FIUX4 




Fluxg^m 


FluX24^m 




(°, J2000) 


(°, J2000) 


(mjy) 


(mJy) 


(mJy) 


(mJy) 


(mj) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


(mJy) 


K194 


22.33638 


-73.557,54 


0.0014 


0.016 


0.47 


0.55 


0.43 


0.27 


0.18 


0.10 


0.047 




K210 


22.42660 


-73.57045 


0.13 


0.35 


0.72 


1.0 


0.66 


0.34 


0.20 


0.13 


0.095 




S213 


22.37980 


-73.56195 


4.3 


2.6 


1.2 


0.79 


0.51 


0.45 


0.30 


0.21 


0.340= 




K225 


22.38890 


-73.56236 


6.4 


3.8 


1.9 


1.1 


0.66 


0.29 


0.16 


0.11 


0.1605 




K232 


22.38700 


-73.56084 


2.0 


1.3 


0.64 


0.37 


0.26 


0.13 


0.062 




0.160=5 




S236 


22.36260 


-73.55448 


9.3 


5.5 


2.66 


1.5 


0.94 


0.39 


0.24 


0.15 


0.12 


9.0O5 


K248 


22.38330 


-73.55859 


1.9 


1.2 


0.60 


0.34 


0.22 


0.13 


0.065 




0.035 




K261 


22.25457 


-73.53501 






1.9 


2.30 


1.66 


0.96 


0.64 


0.40 


0.25 




S293 


22.36570 


-73.54987 


2.1 


1.6 


1.0 


0.68 


0.61 


0.32 


0.22 


0.16 


0.12 




K296 


22.32630 


-73.54260 


1.8 


2.4 


2.6 


2.4 


1.6 


0.77 


0.47 


0.32 


0.24 




K340 


22.36068 


-73.54176 


0.24 


0.22 


0.16 


0.11 


0.082 


0.100^5 










K348 


22.30930 


-73.53139 


4.7 


10 


18 


26 


16 


8.0 


4.6 


2.9 


1.8 




S394 


22.36515 


-73.53374 


0.94 


0.58 


0.41 


0.23 


0.15 


0.14 


0.092 




0.044 




S406 


22.36790 


-73.53262 


2.5 


1.9 


1.2 


0.78 


0.57 


0.42 


0.31 


0.21 


0.15 




S411 


22.32367 


-73.52343 








0.016 




0.10 


0.076 




0.051 




K421 


22.31828 


-73.52143 






36 


36 


22 


10 


6.38 


3.8 


2.2 




K441 


22.29871 


-73.51362 






15 


30 


26 


18 


13 


9.6 


7.7 


0.68 


K442 


22.49467 


-73.54764 






0.066 


0.067 


0.12 


0.29 


0.18 


0.13 


0.073 




K444 


22.34401 


-73.52120 






0.058 


0.11 


0.096 


0.16 


0.13 


0.078 


0.073 




S466 


22.47845 


-73.54208 






19 


26 


17 


8.5 


4.7 


3.3 


2.0 


1.3 


K481 


22.49347 


-73.53891 






0.026 


0.060 




0.17 


0.10 




0.040 




S486 


22.39100 


-73.51968 






1.0 


1.0 


0.64 


0.30 


0.19 


0.14 


0.040 




K701 


22.39838 


-73.51712 






10 


12 


8.6 


3.8 


2.32 


1.6 


0.88 




G109 


22.36616 


-73.57540 






0.016 






0.070 


0.052 




0.18 




G133 


22.28723 


-73.65713 


0.0080 


0.024 


0.2S^° 


0.5290 


0.6300 


0.53 


0.37 


0.18 


0.18 




G150 


22.42516 


-73.57831 


0.0003 


0.0026 


0.019 






0.0548 


0.0346 




0.0075 




G211 


22.27448 


-73.54372 






0.20 


0.29 


0.27 


0.29 


0.22 


0.11 


0.21 




G262 


22.33369 


-73.54814 


0.0016 


0.0093 


0.014 


0.029 




0.073 


0.082 




0.17 




G275 


22.40597 


-73.55920 


0.0033 


0.013 


0.025 


0.054 


0.053 


0.15 


0.15 




0.46 




G343 


22.34403 


-73.63817 


0.0073 


0.018 


0.022 


0.036 




0.028 


0.031 




0.073 




G372 


22.39943 


-73.54250 


0.014 


0.036 


0.05690 


0.062^0 


0.06500 


0.22 


0.16 




0.66 




U364 


22.28734 


-73.52486 








0.044 


0.27 


0.057 


0.052 




0.072 


0.52 


U703 


22.29610 


-73.51102 






0.87 


0.65 


0.52 


0.34 


0.29 


0.22 


0.20 


0.63 



Flux considered upper limit with 90% certainty 
Flux considered upper limit with 95% certainty 



Table 5. Young Stellar Object Parameters. For those sources with optical counterparts, we list the number of 
visually identifiable optical sources. Where possible, we determine masses of the optical sources; others either have no 
good optical photometry or are so young or of such low mass (M < 0.5 Mq) that their masses cannot be determined 

from our evolutionary tracks (See Figure [3] (b)). Mopt is the sum of the optical masses which we can determine. 

Sources with Stage listed as U are unclassified probable YSOs. 



Name 


G07^ 








Mopt 






j^ave 












Evol. 








(cpd'^ ) 




(Mo) 


(Mq) 


(Mq) 




(lO^Lo) 


(lO-^Me/yr) 


(lO-^Mo/j/r) 


(Mo) 


(Mo) 


Stage 


Y090 'f' 


2 


> 10' 


5.39 


2 


0.5-i 


7.56 


0.23 


1.223 


0.496 


20.56 


8.535 


0.1407 


0.0637 




Y096S 


51 


> 10' 


0.31 


69 


ij 


8.04 


0.01 


2.475 


0.013 


1.271 


0.035 


0.2276 


0.0091 




Y118 




1' 


0.22 


6 


0.8 


7.37 


0.09 


0.729 


0.090 


33.770 


3.291 


0.7915 


0.1824 




Y142'^ 


36 


1 


0.02 


120 




7.40 


0.03 


2.016 


0.023 


14.830 


0.212 


0.2022 


0.0040 




Y143'' 




1 


0.39 


1 




7.45 


0.00 


2.257 


0.000 


2.214 


0.000 


0.1467 


0.0000 




Y148 ' 







1.19 


9 




8.17 


0.08 


2.205 


0.010 


3.711 


0.175 


0.7449 


0.1132 




Y149 




N/A' 


0.03 


37 




4.19 


0.06 


0.442 


0.017 


3.572 


0.193 


0.1828 


0.0087 




Y162 


53 


3 


0.73 


15 




6.67 


0.18 


2.581 


0.219 


0.147 


0.025 


0.3671 


0.0425 


l/ll 


Y163 




N/A' 


0.04 


128 




3.64 


0.01 


0.143 


0.008 


3.280 


0.190 


0.1486 


0.0019 


l/ll 


Y170S 


46 


4' 


0.39 


6 




8.23 


0.09 


2.541 


0.166 


4.744 


0.325 


0.1472 


0.0172 




Y171 




7' 


0.43 


35 


0.4-' 


10.95 


0.11 


11.840 


0.353 


33.170 


1.948 


0.2325 


0.0148 




Y174'i' 


38 





0.01 


263 




7.63 


0.01 


2.439 


0.009 


9.728 


0.095 


0.4575 


0.0039 




Y179'' 







0.92 


1 




7.45 


0.00 


2.257 


0.000 


2.214 


0.000 


0.1467 


0.0000 




Y19d ' 




2' 


2.09 


9 


0.6j 


7.74 


0.07 


2.741 


0.070 


0.829 


0.132 


0.2822 


0.0246 


l/ll 


Y197'' 







0.87 


3 




8.11 


0.16 


2.556 


0.137 


2.346 


0.144 


0.5644 


0.1899 


■ 


Y198 




2' 


0.59 


5 




6.85 


0.09 


0.409 


0.018 


18.702 


5.690 


0.7289 


0.1602 




Y206S.'= 


49 


2' 


0.90 


11 


1.4j 


25.87 


0.38 


23.680 


0.354 


103.900 


7.538 


0.5077 


0.0098 




Y223'' 




> 10 


0.24 


18 




6.37 


0.02 


1.202 


0.015 


0.002 


0.003 


0.1762 


0.0110 




Y237S.'= 




2 


0.49 


42 




5.80 


0.04 


1.422 


0.074 


0.110 


0.009 


0.1266 


0.0052 


l/ll 


Y240f 




3 


5.38 


2 


I' 


8.78 


0.38 


4.037 


0.609 


0.533 


0.137 


0.0031 


0.0013 




Y25it,'i,fc 


1 


l' 


4.38 


4 




11.12 


0.90 


14.550 


4.471 


64.900 


23.480 


0.1505 


0.0254 




Y255 




2' 


0.77 


27 


3.55j 


6.50 


0.03 


1.316 


0.013 


< 0.001 


< 0.001 


< 0.0001 


< 0.0001 


ll/lll 


Y264 




7' 


0.32 


10 




6.57 


0.00 


1.295 


< 0.001 


0.000 


0.000 


< 0.0001 


0.0000 




Y271 




3' 


0.05 


31 


2 


4.58 


0.02 


0.367 


0.005 


0.065 


0.012 


0.0001 


< 0.0001 




Y283 


52 


1 


0.63 


20 


3 


4.51 


0.02 


0.349 


0.006 


0.001 


< 0.001 


0.1072 


0.0043 




Y285S''''* 


34 


7' 


1.76 


12 




8.56 


0.06 


3.837 


0.079 


1.859 


0.165 


0.4557 


0.0511 




Y288S 




6 


0.23 


8 


1.8j 


11.25 


0.06 


3.516 


0.086 


6.695 


0.253 


0.3750 


0.0367 




Y290 




2' 


0.05 


46 


5.5j 


6.75 


0.01 


1.494 


0.012 


< 0.001 


< 0.001 


< 0.0001 


< 0.0001 


ll/lll 


Y312S.''.'= 


24 


8' 


2.35 


11 


0.65j 


12.41 


0.36 


6.367 


0.429 


68.470 


7.649 


0.4087 


0.0733 




Y326S.'' 


22 


> 10* 


1.46 


12 


3.85j 


8.01 


0.08 


1.548 


0.086 


2.534 


0.114 


0.7441 


0.0793 




Y327S 


54/57 


> 10' 


0.14 


67 


8.2-i 


8.78 


0.02 


4.175 


0.028 


2.969 


0.143 


0.5354 


0.0130 




Y358 




3' 


0.07 


38 


1.2j 


2.58 


0.01 


0.042 


< 0.001 


0.001 


< 0.001 


0.0408 


0.0017 




Y396 




4 


1.49 


6 




2.69 


0.00 


0.056 


0.000 


0.000 


0.000 


0.0159 


0.0000 





Table 5 — Continued 





G07^ 


* opt 


A. 




opt 










env 


rr A/I 






Evol. 






















('10~*Mri /lir) 










Y493 




N/A' 


0.07 


248 




AAA 


0.00 


0.283 


0.001 


3.914 


0.034 


0.1174 


0.0015 


I 


Y700 




N/Ai 


0.24 


14 




4.31 


0.07 


0.205 


0.003 


0.737 


0.169 


0.0449 


0.0043 


I 


Y217'^ 


26 





0.03 


28 




5.68 


0.01 


0.455 


0.008 


0.514 


0.026 


0.2646 


0.0196 


I/II 


Y217eh 







2.66 


2 




9.74 


0.15 


6.020 


0.427 


0.000 


0.000 


0.6202 


0.3101 


Il/lII 




50 


3 


0.74 


17 




6.86 


0.03 


1.260 


0.015 


4.161 


0.696 


0.4697 


0.0748 


I/II 


Y227eh 




3 


4.32 


6 


2j 


8.87 


0.32 


4.894 


0.436 


0.001 


0.001 


0.3339 


0.0997 


II 


Y287'' 


45 


3 


0.07 


43 


2j 


5.36 


0.01 


0.225 


0.002 


4.895 


0.294 


0.2185 


0.0060 


I 


Y287e'^ 




3 


2.59 


1 


2.2j 


12.04 


0.00 


11.550 


0.000 


0.000 


0.000 


0.4850 


0.0000 


II 


Y270i 


29 


> 10 


1.05 


22 


15J 


9.27 


0.14 


7.400 


0.772 


23.900 


4.978 


0.0025 


0.0005 


I/II 


Y270ie'' 




> 10 


2.90 


10 


15j 


8.17 


0.11 


1.207 


0.051 


41.680 


4.025 


0.4940 


0.0498 


I 


Y270o 




> 10 


1.79 


18 


15j 


5.64 


0.01 


0.472 


0.004 


0.174 


0.006 


0.0666 


0.0042 


I 


A270 




> 10 


3.36 


8 


15j 


8.45 


0.00 


3.463 


< 0.001 


0.000 


0.000 


< 0.0001 


< 0.0001 


II 


A270ch 




> 10 


4.95 


9 


15j 


8.45 


0.00 


3.463 


< 0.001 


0.000 


0.000 


< 0.0001 


< 0.0001 


II 


Y340 


5 


2 


0.62 


29 




4.72 


0.04 


0.409 


0.012 


0.739 


0.049 


0.04732 


0.0048 


u 


Y340eh 




2 


2.84 


5 




7.72 


0.00 


2.459 


0.000 


0.000 


0.000 


< 0.0001 


0.0000 


III 


K340 




2 


0.03 


4138 


4j 


















Stellar 


A340 




2 


3.51 


6 


4j 


7.72 


0.00 


2.459 


0.000 


0.000 


0.000 


< 0.0001 


< 0.0001 


III 


A340eh 




2 


3.61 


6 


4j 


7.72 


0.00 


2.459 


0.000 


0.000 


0.000 


< 0.0001 


< 0.0001 


III 



''ID from lGouIiermis et al.l ll2007h . 

^Nopt is the approximate number of optical sources within the 1'.'6 IRAC aperture radius. 
"xliin (cpd) per datapoint 

'^iT'fits is the number of fits with cpd between cpdmin and cpdmin + 0.5 
'^(T values represent the standard deviation of the mean. 

'These sources are likely YSOs based on CMD examination and physical location, but they are not considered well-fit. 

^These sources are fit with optical fluxes set as upper limits. 

''Fit using the error bar method of correcting for PAH emission. 

'Corresponds to one or more PMS star. 

■■These optical masses neglect one or more optical source. 



■"These are in environments of illuminated dust, mostly along ridges or in the tips of "pillars." 
'Outside optical FOV. 
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